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ABSTRACT 



We introduce a new 1-matrix model with arbitrary potential and the 
matrix-valued background field. Its partition function is a r-function of KP- 
hierarchy, subjected to a kind of £_i-constraint. Moreover, partition function 
behaves smoothly in the limit of infinitely large matrices. If the potential 
is equal to X^^^, this partition function becomes a r-function of i^-reduced 
KP-hierarchy, obeying a set of Wi<-algebra constraints identical to those con- 
jectured in |I| for double-scaling continuum limit of {K — l)-matrix model. In 
the case of = 2 the statement reduces to the early established [0 relation 
between Kontsevich model and the ordinary 2d quantum gravity . Kontse- 
vich model with generic potential may be considered as interpolation between 
all the models of 2d quantum gravity with c < 1 preserving the property of 
integrability and the analogue of string equation. 
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1 Motivation and results 



1.1 Introduction 

Matrix models play an increasingly important role in the theory of quantum gravity, be- 
ing a reasonably simple realization of Regge calculus, adequate at least in essentially to 
topological theories. Their conventional application to string theory is based on identi- 
fication of the critical points of peculiar "double-scaling" limit 0, ^, ^ of Hermitean 
multimatrix models [0 with the c < 1 minimal conformal field theories P] coupled to 2d 
gravity. Such models are often supposed |^ to exhaust the entire set of consistent bosonic 
2-dimensional string models, while bosonic strings beyond 2 space-time dimensions are 
presumably unstable . Technical methods, based on the application of matrix models. 



appeared rather successful not only in the framework of perturbation theory, but allow 
one to estimate the entire sums of perturbation series, thus rising the study of string 
models to a qualitatively new level. 

On the other hand, matrix models appeared to possess unexpectedly deep internal 
mathematical structure 0. In fact, if matrix models are used to describe interpolation 
between different critical points (i.e. particular c < 1 string models), the interpolating 
flows are mutually commuting, i.e. the entire pattern of flows constitutes an integrable 
hierarchy (some reduction of Kadomtsev-Petviashvili (KP) system). This integrable struc- 
ture indicates that the somewhat artificial description in terms of matrix integrals should 
possess a more invariant algebro-geometrical description, which could be used for formula- 
tion of the yet unknown dynamical principle of string theory, unifying all available string 
models into distinguished and (hopefully) physically important entity. 

The main obstacles on the way from matrix models to such general principle are: 

(i) the sophisticated intermediate step: the double-scaling limit, 

(ii) the somewhat obscure origin of integrability: it is not a priori obvious, why matrix 
integrals with generic "potentials" are r-functions of (reduced) KP-hierarchies, 

(iii) the non-universal description of different c < 1 string models: some of them 
are nicely unified as different critical points of particular multimatrix model, but there 
is no nice interpolation between critical points of models with the different number of 
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matrices. The naive way out is to use oo-matrix model, but it is no longer described in 
terms of /inzie-dimensional matrix integrals, and, more important, reductions to the case 
of particular multimatrix model is too much singular. 



In this paper we propose a kind of resolution of all these three problems. Namely, we 
argue that there is a new one-matrix model [|11|] (which we call Generalized Kontsevich 
Model (GKM)), which 

(i) for a specially adjusted potential describes properly the double scaling limit of any 
multimatrix model, 

(ii) has the partition function, which is a KP r-function, properly reducible at the 
points, associated with multimatrix models, and satisfies an additional equation, which 
reduces to conventional £_i-constraints for multimatrix models, 

(iii) allows continuous deformation of potential, changing it from the form, associated 
with a given multimatrix model to those corresponding to the others. 

This makes the study of GKM a natural step in the development of string theory. 

In the remaining part of this section we describe GKM and its relation to multimatrix 
models more explicitly. The proofs are presented in sections 2-4. For a more condensed 
presentation of our results see |]lT] . 



1.2 Generalized Kontsevich model 

This paper is devoted to the study of 1-matrix model in external matrix field A, essentially 
defined by the integral over N x N matrix X : 

J'PlA] ^ jdX e-^'-i^W+^-^^ (1.1) 

with arbitrary "potential" V{X) and dX = Y[i^j=idXij. Such integrals arise in considera- 
tion of various problems [0, |T^, |^, [T^, but this time our purpose is to use (1.1) in order 
to define a universal matrix model, which could describe interpolation between all the 
multimatrix models and their critical points, as explained in sect. 1.1 above. Namely, we 
present strong arguments that partition function of such universal model may be identified 
with ZP[M], where 
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with 

and 



U{M, Y) = Tr[V{M + Y)- V(M) - V{M)Y] 



(1.2) 
(1.3) 



U2{M, Y) = lim -U{M, eY), (1.4) 

{i.e. U2{M,Y) is the F^.^ontribution to U{M,Y)). 
After the shift of integration variable 

X = M + Y (1.5) 
the numerator in (1.2) may be rewritten in terms of (1.1): 

|e-^(*^'^)tiy = e^'-[^(*^)-^^^'(*^)l^i^^[V'(M)]. (1.6) 

The denominator in (1.2) is a simple Gaussian integral — a peculiar matrix generalization 
of [V"(/i)]-i/2. 

We shall call zj^-'^[M], defined by eq.(1.2), the partition function of generalized Kont- 
sevich model (GKM). The reason for this is that for the special choice of potential^: 

V(X) = V2(X) = XV3, (1.7) 
eq.(1.2) becomes the partition function of original Kontsevich model [|l^: 



r rIV 0-1/3 TrY^-TrMY^ 



JdY e 



-TrMY^ 



_ -(2/3)TrA/^ J e /.ON 

^dXe-T-^^' ■ ^ ' 

Eq.(1.8) has a natural generalization, which is also a particular case of (1.2): if 

^ We do not specify integration contour in matrix integrals (1.1), (1-2), since all what we are going to 
discuss does not depend on it. To make contact with the literature, note that for potential Vk = ^k+i 
integrals should be over Hermitean or antiHermitean matrices if X + 1 is even or odd respectively. 
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V(X) = Vk{X) = (1.9) 



eq.(1.2) becomes 



We shall argue that Z^^M], defined by (1.10), is identical to the square root of partition 
function of {K — l)-matrix model in the double-scaling limit with 

T„ = -Tr M-", n ^ mod K , (1.11) 

n 

playing the role of generalized Kazakov's time- variables |16| : 



^i^^[M] = \/rif-^^(T„), (1.12) 

where M 



rjf > = hm n / dM,e^ *„™+T.M,M.+i ^^_^3) 

double scaling J 

(the term M^M^+i in the exponent should be omitted; the change from discrete time- 
variables {tn} to continuum {T„} is implicit in performing the double-scaling limit, see 

miiTi). 



Since potential V{X) can be continuously deformed from one K in (1.9) to another, the 
GKM (1.2) as a corollary of (1.12) may be considered as continuous interpolation between 
all the models of 2d gravity with c < 1, described by particular multicritical points of 
particular multimatrix models (let us remind, that for K fixed the p-th multicritical point 
is defined by the condition [|l| that 



all Tn = 0, except for Ti and Tx+p, (1-14) 

and according to (1.11) this is a constraint on the form of the matrix M). 

This is the sense in which the GKM (1.2) may be considered as the universal partition 
function of bosonic string models. We shall not discuss the implications of this fact 
here, instead we shall concentrate on the arguments in favour of this conclusion, which 
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are entirely based on the identity (1.12) (and also motivated by the nice features of 
interpolating function Zj^J'[M] itself — it is usually a KP r-function, subjected to 
constraint, see below). Our strategy is to try to understand as much as possible about 
Z^>[M], which stands at the Lh.s. of (1.12) and compare these results with what is 
conjectured in |1[] about at the r.h.s. 

1.3 Properties of Z^^^ 

In order to motivate our considerations let us remind first, what is known about Z^^^[M] 
in the particular case of i^' = 2. 



Expression (1.2) has been derived in [|T5| as a representation of the generating func- 



tional of intersection numbers of the stable cohomology classes on the universal module 
space, i.e. it is defined to be a partition function of Witten's 2d topological gravity jTSj . 
In (see also [|l^, [19| for alternative derivations) it was shown that as ^ oo 



considered as a function of time variables (1.11), satisfies the set of Virasoro constraints 



4^=^>Zi^> = 0, n>-l, (1.15) 



k odd ^ a + b=2n 

a,b odd and>0 

+ j aTM + ^6^,o-d/dn+2n. (1.16) 

a,b odd and>0 

(Notations are slightly changed as compared to p|. This derivation is partly reproduced 
in sect. 4. 2 below.). Sometimes it may be convenient to consider the appearance of the 
last item at the r.h.s. of (1.16) as result of additional shift of time- variables, 

T„ -> fi^'> = T^- -6n,K+U (1-17) 
n 

however, these T-times are defined in a -ft'-dependent way and do not seem to have too 



^{1} 

ds ' 



much sense. Constraints (1.15) are exactly the equations [|^, ^ |TB[, imposed on yF 
Moreover, it is very plausible that they possess a unique solution, so that (1.15) in fact 
implies that 
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Alternative solution to the constraints (1.15) is given |]T], ^ by a Galilean- invariant KdV 
r-function , and, relying upon the same belief that the solution is unique, one concludes 
that Z}^^ is a r-function of KdV-hierarchy, subjected to additional £i^/-constraint: 

= = r^'\ (1.19) 

/:12}^{2} ^ Q (^ 20) 

Notation r^'^^ reflects the fact that any rxdv may be considered as a KP r-function, 
evaluated at specific points of Grassmannian (see sect. 2. 5 below for details), corresponding 
to the 2-reduction of KP hierarchy. Generically, the r-function of /^-reduced hierarchy, 
r'f^^(T„), is the KP r-function at peculiar points of Grassmannian and possesses the 
following property: 



' (Ti...Tk-1,Tk,Tk+1--.T2K-1,T2K,T2K+1, ■■■) — 

= e(Sn""^^"^VW(Ti...TK-i,0,Ti,+i...T2K-i,0,T2i,+i,...). (1.21) 

This means that the entire dependence of generic r^^^ of all variables TnK (^.e- Tn with 
n = mod K) is exhausted by the simple exponent 

expC^anKTnK) (1.22) 

n 

with time- independent parameters anx- Moreover, this factor is actually absent in the 
case of GKM partition function: a^xiZ^^^] = 0. 

The significance of the just reported results concerning = 2 is two-fold. First, 
they establish the relation (1.7) between two different models of 2d gravity: Witten's 
topological gravity [|T^ , represented by [|I^ , and ordinary quantum 2d gravity 0, ^, ^ , 
represented by T^J . This was the main emphasize of ref.0]. The second implication, and 
it is the one of interest for us in this paper, is the possibility to describe the sophisticated 
double-scaling limit of matrix model in terms of a completely different matrix model (1.2). 
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It is important that is in fact a "smooth" function of as ^ oo, and all the 

properties of are just the same as at finite A^'s, — in variance with the double- scaling 
limit for ordinary matrix models, the limit cxo for GKM is non-singular. Also, to 

avoid misunderstanding, we should emphasize that what we mean under double- scaling 
limit here is not the limit, describing a particular fixed point, but the entire pattern of all 
the flows between various fixed points of a given {K — l)-matrix model, so that partition 
function T^f^ is indeed a function of all time- variables (1-3). 

1.4 The case of arbitrary K 

In order to generalize the identity (1.15) to the case of arbitrary K, i.e. to prove (1.12), 
one can try to prove either an analogue of (1.15) or that of (1.19), (1.20) and compare the 
result with the what is expected for \fr^~^. This will be our strategy below. Namely, 
in sect. 2 a simple and promising formalism is developed, based on interpretation of (1.11) 
as Miwa transformation, and we use it to prove that any Z^^[M], defined by (1.2) and 
(1.11) {i.e. with arbitrary V(X)) is in fact a KP r-function — in particular, satisfies 
bilinear Hirota difference equations. Moreover, in the Grassmannian parameterization of 
KP r-functions the point of Grassmannian is specified by the choice of potential V{X), 
so that multimatrix models with different K are (due to (1.12)) associated with different 
points of Grassmannian. (In terms of the universal module space these points are in fact 
associated with certain infinite-genus hyperellyptic surface for K = 2, and with certain 
infinite genus abelian coverings of degrees K in the general case. In this sense interpolation 
with the changing V(X) is some flow at infinity of the universal module space.) 

What is specific for a given K is that whenever V{X) is a homogeneous polynomial 
of degree K + 1, i.e. if it is of the form (1.9), the KP r-function becomes independent 
of all Tku, i-G- acquires the property (1.21) and may be considered as r-function r^^^ of 
i^'-reduced hierarchy. This will be also explained in more details in sect. 2 below. 

Though any Z^J'[M] defined by (1.2) is a KP r-function, the inverse statement is not 
correct. Z^J'[M] is subjected to the infinite set of constraints, essentially implied by the 
Ward identities for the integral (1.1) ||12|, 
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Tre(A) 



J^lr^^ = (1.23) 



with any e(A). These identities result from invariance of (1.1) under any shift of integra- 
tion variable X ^ X + e(A). Eq.(1.23) seems to provide a complete set of differential 
equations, unambiguously defining (up to a A and A^-independent constant). This 
statement is the actual basis for all of the arguments, which make use of the uniqueness 
of solutions to any constraints. 

As shown in in the particular case oi K = 2 the entire set of such identities with all 
possible e(A) is equivalent to (1.15) and thus implies (1.19) and (1.20). In fact in sect. 2 
we follow a somewhat different logic and give an alternative direct derivation of relations 

Z^:^^ = r, ZW = rW (1.24) 

for any V{X) and so it remains only to prove the analogue of (1.20). Such direct proof 
is presented in sect. 3, moreover it is valid not only for potentials of the form (1.9), but 
for any V{X) : 



\4l ^V"(M) 



d 

+ 



;i.25) 



where /Zj are eigenvalues of the matrix M. 

For potentials (1.9) acquires conventional form of C^^]^ (see eq.(1.28) below) and 
may be easily supplemented by other generators of Virasoro algebra. 

Despite of eqs.(1.24) and (1.25) provide a complete description of partition function of 
GKM, the study of direct consequences of Ward identities (1.23) also seems interesting. 
In particular, it should lead to a direct derivation of the analogue of Virasoro constraints 
(1.15). It also deserves noting that the meaning of constraints on Z^^[M] is that V{X) 
parametrizes only some restricted subset of the universal module space. Only this subset 
is associated with matrix models and interpolations between them. This is also clear 
from the fact that V{X) = J2 V„X" is naturally parameterized by an infinite vector 
{V„,} rather than by a matrix Amn, what would be natural for the entire Grassmannian. 
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Remarkably enough the "dimension" of this subset is just the same as that of the space 
of time-variables (the space of potentials is parameterized by a vector {V„}, while that 
of time- variables — by {T„} or by eigenvalue vector {/in} of the matrix M). This clearly 
suggests that a formalism should exist for so restricted set of r-functions, which would 
treat potentials and times in a more symmetrical fashion (see also discussion in sect. 3. 4)^. 

Implications of Ward identities (1.23) for Z^^[M] with K > 2 are discussed to some 
extent in sect. 4. This is a straightforward calculation, though much more tedious than 
in the simplest case of = 2. We argue that the proper generalization of (1.15) for the 
case of = 3 is given by the following set of constraints: 



W^Zj^^ = 0,n> -2; 

^k>l a+fe=3n '^^ 3a+2 J 

Y: {3k-2)f,,_Mt3n+ E WgLa Ui?^ = 0, a, & > 0, n > -3. 

fc>l+<5„+3,o a+b=3n 3a+l J 

(1.26) 

Here stands for the i^n-th harmonics of the j-th generator of Zamolodchikov's Wk- 
algebra, and the proper notation would be W^^^^^ or W^]^^\ The way to express these 
quantities through free fields is described in [jl|; for example, 



= 3 y kfkifi—-^ — + 3ykfk y y ^ ^ + 

k,l>l 01k+l+3n fc>i a+b=/fc+3n '^^aCJb ^^^^^^^^^ Ol aOl bOl c 

+ y^ aTabTbcTc, a,b,c > 0, a, b,c,k,l ^ mod 3 

a+b+c=—3n 

(1.27) 

and 



^ ' 2ir 

a,b>0 



^Let us also remark that such symmetry exposes itself in the fact that the degree of potential K and 
the order of the multicriticity p + K (fixed by the choice of times) give the theory of minimal matter with 

p-K 
pK 



c = 1 — coupled to 2d gravity, and, thus, p and K should enter the theory on equal footing. 
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+± J2 ,tM^(Ji^}}iii±}ls„,. (1.28) 
2A'.^.-e^, 24A' 

a,b>0 

(Note that in order to make formulas a bit more compact, we expressed all these operators 
through T^^^ rather than T„-variables. These are defined in (1.17), d/dT = d/dT.) The 
analogue of (1.26) for > 3 is a more or less obvious generalization of (1.26), involving all 
the generators wj^\ ^\ wj^^ of WA-algebra with the "coefficients" made from 

T's and {d/dT)'s. The form of these operators is a bit reminiscent of the W-operators, 
introduced in p^ , |22| as ingredients of Ward identities in discrete multimatrix models. 

Eqs.(1.26) (and their analogues for > 3) are obviously resolved by any solution of 
the simpler set of equations, 

WS^W = 0, A; = 2, 3, n > 1 - (1.29) 

(eq.(1.20) is the particular case: k = 2, n = —1). On the other hand (1.26) are equivalent 
to (1.23) for the particular form of potential (1.9), and thus are supposed to possess a 
single solution. Therefore from this point of view it also looks plausible that (1.26) (and 
its counterpart for any K) is simply equivalent to (1.29). Unfortunately we do not provide 
an explicit proof of this equivalence (though we could rely, of course, upon the indirect 
proof from ref.[Q, deducing (1.29) from eqs. (1.19), (1.20), which are proved in sect. 2,3 
without any reference to (1.23)). 

If one believes in this transition from (1.26) to (1.29), we can finally return to multi- 



matrix models. Namely, in [jl| it was suggested that yT^f 

(i) is a r-function of /^-reduced hierarchy, 

(ii) satisfies the constraint 



£iTVrif "'^ = 0, n > -1, (1.30) 
and, corollary of (i) and (ii), 

(iii) satisfies the entire set of eqs. (1.29): 



wSvrif "'^ = 0,k = 2,3,...,K;n>l-k. (1.31) 

Now it remains to use the above-mentioned assertion that this system has a unique solution 
in order to arrive to our identification 
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Z{K} 

and thus to the central idea of this paper 



: universality oi GKM (1.2). 



(1.32) 



1.5 Comments 

Before we proceed to actual derivation of all these statements, let us comment briefly 
on the very status of the constraints (1.29) imposed on \Jv\^ . In order to honestly 
derive such constraints on the lines of |jT6| it is first necessary to find out their counterparts 
(loop equations) in the discrete multimatrix models (these are expressed in terms of the so 
called W operators in 0, |2^), then carefully perform the appropriate reductions (including 
identification of all the K — 1 a priori different potentials), and take the double- scaling 
limit. The actual form of the W-constraints is very similar to (1.26) and therefore it may 
happen that the final constraints would arise just in the form of (1.26) rather than (1.29). 
This could make the relation (1.26) even more important. 

After this general description of our reasoning and results it remains to prove that 

(i) Z}^^ as defined by (1.2) (the limit — > oo is smooth) is indeed a KP r-function; 

(ii) for any V{X) it is subjected to additional constraint (1.25); 

(iii) if V(X) = const ■ X^^^, then in fact becomes a r-function t^^^ of the 
A'-reduced KP-hierarchy; 

(iv) in this case also the entire set of constraints like (1.15) for A' = 2 and (1.26) for 
K = 3 are imposed on Z^\ 

The proof of the points (i),(iii) is given in sect. 2, (ii) is discussed in sect. 3 and (iv) 
— in sect. 4. While (i),(ii),(iii) are proved in full generality — for any potential V(X), — 
the universal form of the constraints (iv) is not completely understood. Moreover, the 

XK+l 

proof of (iv) for the case of V(X) = , as presented in sect. 4, is not complete even 

K -\- 1 

in the simplest case of i^' = 3. Also the proof of (iv) is not independent of (iii), since 
eq.(1.21) is used. The deep understanding of these Woo-constraints and their universal 
formulation remains the first obvious problem for GKM to be resolved in the future. The 
second obvious problem is to find out the algebro-geometrical origin of the model (1.2) 
itself. 
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2 Partition function of the GKM as KP r-function 



2.1 Short summary 

The purpose of this section is to prove that: 

(A) The partition function zj^\M], if considered as a function of time- variables 

T„ = -trM-- , (2.1) 
n 

is a KP r-function for any value of N and any potential V[X]. 

(B) As soon as V[X] is homogeneous polynomial of degree K + 1, zjf\M] is in fact 

a r-function of K-reduced KP hierarchy. Moreover, actually, — — — = 0. 

oTnK 

In order to prove these statements, first, in sect. 2.2 we rewrite the r.h.s. of eq.(1.2) 
in terms of determinant formula (2.21) 

Z^'lA/l = '^^''2(tf '■^ = 1.-.'V- (2-2) 
We show in sect. 2.3 that any KP r-function in the Miwa parameterization does have the 
same determinant form. Finally, as a check of self-consistency, we prove in sect. 2.4 that 
any determinant formula (2.2) with any set of functions {(j)i{fi)} satisfies Hirota difference 
equation. 

A brief description of Grassmannian (or Universal module space) language is given 
in sect. 2.5. The point of Grassmannian corresponding to r-function (2.2) is defined by 
the infinite- vector function {0j(/i)}. Remarkably, the results of ref. which concern 
identification of the points in Grassmannian, associated with multi-matrix models, are 
immediately reproduced from the explicit expressions for the basis {(j)i{fi)}. This exhausts 
our discussion of (A). 

The main thing which distinguishes matrix models from the point of view of solutions 
to the KP-hierarchy is that the set of functions {0i(/u)} in (2.2) is not arbitrary. Moreover, 
this whole infinite set of functions is expressed in terms of a single potential V[X] {i.e. 
instead of arbitrary matrix Aij in (j)i{fi) = J^^ijf^'' we have here only a vector Vi or 
"^[fA = Vi/i*). This is the origin of £_i and other W- constraints (which in the context 
of KP-hierarchy may be considered as implications of i2_i). All these constraints are 
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in fact contained in the Ward identity (1.23). While the detailed discussion of these 
constraints is postponed to the sect. 3,4, in subsect.2.6 we discuss another implication of 
(1.23) arising in the case of polynomial V[X]. Namely, whenever V[X] is a homogeneous 
polynomial of degree K + 1, 0j+x(/i) in (2.2) can be substituted by and in such 

points of Grassmannian KP r- function is known to possess the property (1.21) and thus 
may be considered as a r-function t^^^ of fC-reduced KP hierarchy. In sect. 2. 7 we prove 
that the factor (1.22) is actually absent in GKM. This is all to be said below about (B). 

2.2 From GKM to determinant formula 

Numerator of (1-2). We begin from evaluation of the integral (1.1): 



If eigenvalues of X and A are denoted by {xi} and {Aj} respectively, this integral can be 
rewritten as 



V/v stands for the volume of unitary group U{N) and is unessential in what follows; 
A(X) and A(A) are Van-der-Monde determinants, e.g. A(X) = Y\iyj{xi — Xj). The 
transformation, leading from (2.3) to (2.4), is a trick, familiar from the study of multi- 



dwell upon this point. For what follows it is very important that the term TrAX in 
(2.3) is linear in X (be it instead, say, TrAX^ , the polynomial factor A(X) at the r.h.s. 
of eq.(2.4) would be substituted by A^{X)/A{Xp) and our reasoning below would be 
unapplicable) . 

The r.h.s. of (2.4) is proportional to the determinant oi N x N matrix 




(2.3) 




(2.4) 



matrix models in the formalism of orthogonal polynomials |^, and we do not 



A-i 



(A)det(,j)Fi(Aj) 



(2.5) 



with 
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F,+i(A) = ldx x^e-^(^)+^^ = (|^)^^i(A). (2.6) 

Note that 

F,(A)=^S[A] . (2.7) 

If we recall that in GKM 



A = v'(M) (2.8) 
and denote the eigenvalues of M trough {/Xj} , eq.(2.6) acquires the form: 

where 



!.,(/.) =F,(V'(/.)). (2.10) 
Denominator. Proceed now to the denominator of (1.2): 

lP[M] = JdX e-^^(^'^). (2.11) 

Making use of [/(A'^)-invariance of Haar measure dX one can easily diagonalize M in 
(2.11). Of course, this does not imply any integration over angular variables and provide 
no factors hke A(X). Then for evaluation of (2.11) it remains to use the obvious rule of 
Gaussian integration, 

j dX e" ^'i^'i^^^ ^ [] j-2.12) 

(an unessential constant factor is omitted), and substitute the explicit expression for 
Uij{M). If potential is represented as a formal series, 

oo 

V{X) = J2 K^"' (2-13) 

n=l 

(and thus is supposed to be analytic in X at X = 0), eq.(1.4) implies that 
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[/2(M,X) = ^(n + l)V„+i ^ TrM'^XM'X}, 

n=0 I a+b=n—l 



and 



n=0 {a+b=n-l J n=0 A^i A*J' 

Coming back to (1.2), we conclude that 



Here = [[/..] V2eV(w)-MiV'(/.i)^ ^_g_ 

= [V"(/i)]'/'e^^''^~''^'^''^ (2.15) 
The product of s-factors at the r.h.s. of (2.16) can be absorbed into ^-functions: 

Z<n[M] = t!^^, (2.16) 

where 

$,(//) = s(/x)|.,(/x). (2.17) 

Kac-Schwarz operator. From eqs.(2.6),(2.10) and (2.15) one can deduce that $»(//) 
can be derived from the basic function $i(//) by the relation 



$,(/x) = [V"(/x)]V2y(^ + ^)-ie-^('^'^)ciy = A-;(^)$,(/x) , (2.18) 
where U{ii,y) is defined by eq.(1.3) and A^yy{ii) is the first-order differential operator 



v"{^l) dfx 2[v"(//)]2 ■ 
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In the particular case (1.9) 



Id k-1 

coincides (up to the scale transformation of /i and As^K}{^J)) with the operator which 
determines the finite dimensional subspace of the Grassmannian in ref . [^] (we shall return 
to this point in sect. 2. 5). We emphisize that the property 

$,+!(/.) = A{v}{^^)U^^) (F,+i(A) = ^i^.(A)) (2.20) 

is exactly the thing which distinguishes partition functions of GKM from expression (2.2) 
for generic r-function, 

rif\M] = M|M^, (2.21) 

with arbitrary sets of functions (f)i{n)- In the next section we demonstrate that the quan- 
tity (2.21) is in fact a KP r-function in Miwa coordinates. Implications of additional 
constraint (2.20) will be discussed in sect. 2. 6 and sect. 3,4. 

N -dependence. Before we proceed to discussion of Miwa coordinates, two more things 
about the formula (2.21) deserve mentioning. First, the entire set {(f)i{fi)} is naturally A^- 
independent and infinite. It is reasonable to require also that </)j's are linear independent 
and 

0,(/i)/0i(/i) = Q/i'-^(l + o(-)), (2.22) 

with /i- independent Cj (this is true for functions (2.7)). Then the set {(/)j(/i)} is identified 
as projective coordinates of a point of Grassmannian, see sect. 2. 5. The r.h.s. of eq.(2.21) 
naturally represents r^'^[M] for an infinitely large matrix M. In order to return to the 
case of finite it is enough to require that all eigenvalues of M , except fii, . . . , fiN, tend 
to infinity. In this sense the function r^'^'-'^[M] in (2.21) is independent of ; the entire 
dependence on comes from the argument M: A^ is the quantity of finite eigenvalues of 
M . As a simple check of consistency, let us additionally carry /ijy to infinity in (2.21), 
then, according to (2.22) 
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detjv0i(Atj) ~ cn4>i{ij^n){iJ'n)^ • detjv-i0i(Aij) • (1 + o{1/ij,n)) 

and 

An{M) - (/.^,)^-lA^^-l(M)(l + o{1/i^n)), 

so that 



rlfi • [c^0i(/.iv)](l + 0(1//.^)). (2.23) 



This is the exact statement about the A'"-dependence of tn- Actually in GKM CAr0i(/x) = 1 
at // — > oo . In what follows we often omit the subscript N. 

Multiplicities. The second remark concerns the form of the r.h.s. of (2.21), when some 
eigenvalues of M coincide (see also eq.(2.30) below). If eigenvalue /Xj appears with the 
multiplicity Pi , eq.(2.21) looks like 

det[0i(/ij)A|v}0i(/ij) ... A^{;]Vi(/^i)] ,^ ^ 

n,>M-,,r- ■ 

Notation in this formula is not very transparent: it is implied that the matrix in the 
numerator has rows of the form 

0j(//l), ^{v}0i(/^l), - ,^{v}Vi()"l), 

MfJ'2), A{V}MfJ'2), ■■■ ,^{^Vj(/^2), MfJ's), ^{V}0i(/^3), ••• , 

where operator A{v} is defined by eq.(2.19). Expression (2.24) will be used in sect. 2. 4 
in the proof of Hirota identity. Note also that if some //j — > oo, this may be treated as 
vanishing of the corresponding multiplicity, pi — 0. Such //j obviously do not contribute to 
(2.24). Thus, the value of N (the finite size of the matrix M in GKM) may be interpreted 
as the number of /x's which appear with non- vanishing multiplicities. 

2.3 KP r-function in Miwa parameterization 

Generic KP r-function is a correlator of a special form: 
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T^{Tn} = (0| : : G\0) (2.25) 

with 

J{z) = ilj{z)i>{z); G = : exp Qmni'mi'n ■ (2.26) 

in the theory of free 2-dimensioiial fermionic fields ip{z), ip{z) with the action / ipdip. 
Vacuum states are defined by conditions 

V'niO) = n< , V'niO) = n > (2.27) 

where ij{z) = Ez V^n-^" dz^^'^ , ip{z) = Ez ^n^"""^ dz^/\ 

The crucial restriction on the form of the correlator, implied by (2.26) is that the 
operator : : G is Gaussian exponential (exponent is quadratic in fields), so that 

insertion of this operator may be considered as modification oi ip — ip propagator, and 
Wick's theorem is usually apphcable. Namely, the correlators 

Cg{M,{>^j}) ^ {0\Yl^{^i^)^{X^)G\Q) (2.28) 

i 

for any relevant G are expressed through the pair correlators of the same form: 

Cg{M, {A,}) = det(,,)CG{/^i, A,}. (2.29) 

Operator G in (2.28) can be safely substituted by the entire : e^^""^" : G , but we do 
not need this for our purposes below. Instead with the help of Miwa transformation we 
shall express : e^^""^" : in (2.25) through insertions of fermionic operators. Then (2.25) 
acquires the form of (2.28) and after that the application (2.29) provides a representation 
for T-function in determinant form, to be compared with eq.(2.21) for partition function 
of GKM. 

There are slightly different forms of Miwa transformation. Usually one writes 

Tu^-EpA (2-30) 



21 



over Riemann sphere with coordinate /i and 

singular function = ^iPiS{ii — /li). We use instead the transformation (2.1): 

T^^-TrM-^^-J2^ , (2.31) 

interpreting /i^s as eigenvalues of the matrix M , all coming with unit multiplicities Pi — 1 
(while Pi > 1 in (2.30) may be interpreted as result of coincidence of Pi eigenvalues). The 
form (2.30) of Miwa transformation is preferable from the point of view of invertibility 
and also is convenient for the proof of Hirota identities (sec sect. 2. 4 below). 

The simplest way to understand what happens to the operator eS^"''" after the sub- 
stitution of (2.31), is to use the free-&oson representation of the current J{z) — d(p{z). 
Then E T^Jn = | -—^^n \ = Yl fil^i)^ 

.gEiV(/^i):= A -TT:e'^('''^ . (2.32) 

In fermionic representation it is better to start from 



instead of (2.31). Then 



: e^^"^" := , , n^(/^0^(/^0 ■ (2-34) 

In order to come back to (2.31) and (2.32) it is necessary to shift all /ij's to infinity. This 
may be expressed by saying that the left vacuum in (2.25) is substituted by 

(A^l ~ (0|^(oo)^'(oo)...V5(^-^)(oo). 
The T-function (2.25) can be represented in various forms: 



T^[M] = (0| : e^^"-^" : GjO) = i^{M)-\N\ JJ : e'^^'^'^ : G\Q) = 

i 

= -lim — ''^-^ — p — (0 j|VG"i)V'G"i)G' 0) = Jim „ ..^ . . . . , 

(2.35) 
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and applying the Wick's theorem (2.29) we obtain: 



-^[M] = hm det,,)^#l4 = (2.36) 

with functions 



0,(//) ~ (0|V;(^-'noo)V'(/^)G|0). (2.37) 

Thus, we proved that KP r-function in Miwa coordinates (2.2) has the determinant form 
(2.1). 

Below we shall need a more detailed expression of the right hand side in the eq.(2.35), 
when all points jli tend to the same value 1/t with fixed t ^ 0. Prom (2.35) one can 
obtain: 



T^m] ^ ^ii^-^^det l-^^dr{om-'m^ij)G\o)\ = 

In the limit of t — > we obtain the formula (2.37) with 

1 (Q\ ■ p^^('"""-*")-^" ■ G\0) 
Un) = lim t) = lim dl^ — ^ . (2.39) 

Eqs.(2.38) and (2.39) will be used in sect. 3. 2 in one of our derivations of the con- 
straint. 

One more remark is in order now. Prom the form of the operator expansion for 
■0(oo)'0(a*) it is clear that the function 0i(/x) , as given by (2.39), behaves as iJ'~^{l-\-o{l/ ji)) 
when oo, in accordance with (2.37), and all the functions 0i(/x) are independent. This 
last statement should be more or less clear from the fact, that the set {0i(/i)} contains all 
the information about fermionic propagator CG(/i,/i) , which is exactly the same as the 
information in operator G , which should be of the form (2.26), and defines the actual 
fermionic action. To make this mutual independence of 0's even more transparent we prove 
in the next section 2.4 that Hirota identities are satisfied by the quantity (2.36) with any 
set of functions {0j(/x)}. The possibility to choose {0i(/x)} in any way is important for 
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us to argue without going into any more details, that partition function of GKM, which 
was expressed in determinant form (2.36) with specific choice of 0j = , is indeed a KP 
r-function. 



2.4 Hirota equations for r-function in Miwa coordinates 

The usual form of bilinear Hirota equations, which are the defining equations of KP 
r-functions, is: 

oo 

^P,(-2|/)P,+ipT)eE.^'^^«V ■ r = 

1=0 

where Dt are Hirota symbols, D = {Dt^, hDx^, ■ ■ •) and Vi are Schur polynomials. 



Note that these equations are in fact more than KP equations themselves, which 

d'^ log T 

describe evolution of the functions Mi(T„) (u(T„) = U2(T„) = — — and all other Ui(T„) 

oil 

d d 

can be determined from the relations — — -— — -logr = with L = d + J2t '^i+i^ * 

oTi oTn 

d 

and d = -— see (3.35) and (3.36)) rather than r-itself. Indeed, any transformation 

oil 

t{T^) - r(T„)e^(^-^-)+^^^ (2.40) 

with dH/dTi = and H = const does not change Ui{Tn), i.e. H and H can not be 
fixed by the entire set of KP equations. They are, however, fixed by Hirota equations (up 
to linear function H + TiH = J2 i>nTni = const). This remark is important for the 
subject of matrix models when ever interpolation between critical points is considered. 
For example, one can use the representation 

r(T„) = exp{ I^^Ti - x)m(x,T2 . . .)dx + H{T2,T^ . . .) + TiH} (2.41) 

and at any given multicritical point the functions if, H are unessential. However, they 
are very important for interpolations, and, in particular, for the relevance of Virasoro and 
W-constraints. 

Our first task in this section is to rewrite the Hirota equations in Miwa coordinates 



(2.1). Since this is a widely known transformation (see, for example |2^), we just cite 
here the answer: Hirota equations in Miwa coordinates state that 
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(A*a - IJ'b)r{Pa, Pb, Pc + l)r(pa +l,Pb + 1, Pc) + 
+ {fJ'b - fJ'c)r{Pa + l,Pb,Pc)T{Pa,Pb + l,Pc + 1) + 
+ (/^c - tJ'a)T{Pa,Pb + l,Pc)T(Pa + l,P6,Pc + 1) = 0, (2.42) 

where r-function is expressed through the Miwa variables (2.30) according to eq.(2.24) 
and three arbitrary eigenvalues from the set {//j} with corresponding multiplicities are 
chosen. 

The second task is to prove that t[M] = A(M)~Met 4>i{nj) with any {0j(/x)} satisfies 
(2.42). 

All we need to derive the desired equations is the well known Jacobi identity for the 
determinants. For any (N + 2) x (N + 2) determinant J we denote by J the 
determinant obtained from J with s rows ii, ... , ig and s columns ji, ■ ■ ■ ,js omitted. 
Then the Jacobi identity reads 

Let us consider TAr({/i}) = TN^fJ^i, • • • , /xat) and divide the given ( a priori different) set of 
eigenvalues /ii, . . . , /in into the L "clusters" of sizes pi, . . . ,pl : 

where J2jPj = N. Then one should introduce two additional eigenvalues /xat+i, ij>n+2 and 
apply the Jacobi identity for i = + 1, j' = A'^ + 2 to the function tn+2{{iJ'}, A^at+i, A*Ar+2)- 
Then simple calculation gives the following system of equations: 

TN+2{{iJ'}, /^AT+l, /^Ar+2)TAr({/^}) = 

= {tn+i{{ij>}, /Xiv+2)TAr+i({M, IJ'N+i) - (2.43) 

- TAr+l({/^}, HN+l)fN+l{{l^}, I^N+2)} 

where f^v+i denotes some new r-function which is obtained from the given tjv+i by a 
change of the last row: (pN+i — ^ (pN+2 ■ Now let all the eigenvalues in each cluster tend 
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to the values jii, . . . , (11 respectively and jJ-N+i l^a, l^b where and /if, belong 

to different clusters. Then eq.(2.43) acquires the form (in the notations as in eq.(2.42)): 



(yUa - fJ^b)r{Pa + l,Pb + 1, Pc)r{Pa, Pb, Pc) = 

(2.44) 

= riPa, Pb + 1, Pc)r{Pa + 1, Pb, Pc) - r{Pa + 1, Pb: Pc)r{Pa: Pb + 1, Pc) , 

where pc describes some arbitrary third cluster different from the previous ones. One can 
multiply this equation by the factor (^"'-^^'-P*^ ) ^^^^^ write down the couple of two 

r{Pa,Pb,Pc) 

another equations obtained by cyclic permutations among indices (a, b, c). The sum of 
these three equations coincides with eq.(2.42). 

Another interesting expression can be derived from the eq.(2.44) as follows. Let us 
make the shift Pa Pa — ^, N ^ N ^ 1 and put pb = (this means that we have the 
cluster with the single element /i(, in the corresponding r- functions in eq.(2.44)). Then in 
the limit /ib — > we have r(pa, Pb + 1) ^ T{Pa + 1) (notations of other clusters will be 
omitted since they don't changed) and equation now takes the form 



TN+l{Pa + l)TN-l{Pa " 1) = 



d d 

= Tiv(Pa) lim T^fj^iPa, fJ'b) " Tiv(Pa) lim ^T^iPa, fJ-b) ■ 



Simple calculation gives further that 



lim -^TN{Pa, IJ-b) = —-^TN{Pa) 



and now we obtain 



d T\X) 1 

PaTN+l{Pa + l)TN-l{Pa " 1) = T%{pa)^ ^ . (2.45) 

This equation holds for every function in the form (2.21). More concrete expression can 
be obtained for the GKM partition function J^j^^ defined by eq.(2.9) in the terms of A's 
variables (see eq.(2.8)). For this quantity the following relation holds: 
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Therefore we have 



Pa^i+'l(Pa + l)^i-'l(Pa-l) 



We should remark that in the case of pi = N, pa = 0, a > 1 (A^ = A) this equation 
reduces to the Toda-chain one. Indeed, now J^n'^Po) has very simple determinant form 

It is just determinant of matrix with specific Toda-chain symmetry of constant entries 
along anti-diagonal From the other point of view, matrix integral (1.1) reduces in 



the present case (of proportional to unit matrix A) to standard discrete model partition 
function which is well-known to correspond to Toda-chain hierarchy |2^, |l^ . If we return 
to the original description in the terms of /x's (see eqs.(2.16), (2.24)), then our final 
equation is 

rN+l{Pa + l)TN-l{Pa-l) ^^2..y ^ (2 46) 

A/" = f n - f^nY-A n ([v"(/i„)]^/^e^('^")-'^"^'('^"))'" . 

\m>n / 71=1 

We shall return to eq.(2.46) in sect. 3. 3 in the context of the string equation. 

2.5 Grassmannian description of KP r-function 

In this section we present the necessary material about infinite dimensional Grassmanni- 
ans and their connection to integrable hierarchies. The details and proofs can be found 



in P, RO, P . 



Definitions. Let us consider the infinite dimensional Hilbert space H realized as the 
space of Laurent series 0(yu) on the unit circle S^{fi G S^) on the complex plane. This 
space is naturally represented as the direct sum H = H^Q)H_ , where H+{H_) is spanned 
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by the vectors {/i"}, n > {{n^"'}, n > 0), n G Z. These vectors form the orthogonal 
basis in H: (/x", /x™) = 6nm- 

The infinite dimensional Grassmannian Gr is the set of subspaces W & H (Grassman- 
nian point) satisfy the two conditions: 

1) the orthogonal projector pr+ : W — > if+ is Predholm operator, i.e. the kernel and 
co-kernel are finite dimensional spaces; 

2) The orthogonal projector pr_ : W — > if_ is compact operator. 

This definition corresponds to Sato's Grassmannian rather than Segal- Wilson one, 
where should be convergent in some neighborhood of oo . 
The index of pr+ is called the relative dimension of W : 

ind pr+ = dim(A;er pr+) — dim{coker pr+). 
Grassmannian is decomposed into the direct sum of the connected components: 

Gr =© Gr„, (2.47) 

n 

where Gr^ consists of VF's with the relative dimension n. For our purposes it is sufficient 
to deal with Gvq. 

It is convenient to represent linear operators acting in H as block matrices with respect 
to the decomposition © H_ : 



a 



, a e GL{oo), (2.48) 



a b 
c d 

where the operators a, b, c, d act as follows: a: if+ b: and so on. 

Clearly, b and c should be compact operators. 

For applications to integrable hierarchies the following commutative subgroup F C 
GL{oo) is of importance. It consists of the mappings from to non-zero complex numbers 
and acts on H by multiplications of the Laurent series. There are two natural commutative 
subgroups r+ and F_ in F : the elements of F_|_(F_) can be analytically continued inside 
(outside) 5*^. Their elements 7+(/i) G F_|_ and 7_(/i) G F_ can be parameterized through 
the infinite sets {T^} and {T^} {k > 1) of independent variables ("times"): 
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7+(/i) = expiYl Tkfi"), 74/") = exp(E Tkf^^')- (2-49) 

k>l k>l 

Let {(pn} > 0) be a basis in W. It is called admissible if the operator transforming /i"^ 
to prj^{(f)n) in has a well-defined determinant. 
We can write 

oo oo 

0„ = + (2.50) 

k=l k=l 

where 0+ has a determinant, 0_ is compact. In particular, for W G Gtq it is always 
possible to choose {(j)^^^)mn = ^mn- In this case we call this basis {0'-'^^} canonicaHoi W. 
Generally, the convenient choice for admissible 0+ (to be used below) is a lower-triangular 
matrix with unit diagonal elements: (0+)nfc = if n < /c. In this case 

0fc(/x)=^'=-i(l + O(/i-i)), A;>1 (2.51) 

(compare to (2.22)). 

The r-function of KP hierarchy tw is defined as determinant of the orthogonal pro- 
jection •y+W : 

rw{{Tk}) = det(pr+ : i+{{n})W ^ H+), (2.52) 

or, more explicitly, 

rw{{Tk}) = det(l + a-^bB), (2.53) 

where 7+ is represented in the block form (2.48), B = 0_(0+)^^ (see (2.50)). 

Given a Riemann surface of finite genus one can associate with it a point in Gr by 



constructing the admissible basis according to the well-known rules In this case 

the formal Laurent series for the basis vectors are in fact convergent. Our situation is 
different: our basis vectors are asymptotic (formal) series which do not need to converge. 
This should be interpreted as adequate generalization to infinite genus Riemann surfaces. 
In this sense we often identify Gr with the universal module space (of line bundles over 
Riemann surfaces with punctures). 
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Instead of W one can work with the equivalent space ■y-{{Tk})W , the r-function being 
changed in the trivial way: 

oo 

T^.wiiTk}) = exp(^ kTkfk)Tw{{n}) . (2.54) 

k=l 

Let us note that in Miwa variables (2.30) the "evolution factor" 7+(/u) (2.49) looks as 
a product of pole factors: 

7+(/^)=n(i--)""^ • 

The various reductions of the KP hierarchy can be described in these terms as follows. 
Let /(/i) be a Laurent series in /x. Those W^s which satisfy the condition 

f{l2)W C W (2.55) 

give rise to solutions of the "/-reduction" of KP. In particular, the choice f{fi) = fi^ 
corresponds to the well known /^-reductions (KdV, Boussinesq, ...) which are in fact 
associated with A^-i Lie algebra [^]. In this case the r-function does not depend on T„ 
with n = mod K (modulo trivial exponential factors linear in times like that in (2.54)). 
An admissible (but non-canonical) basis can be chosen in such a way that 

(j)K+n{^) = fi^(pn{fJ'), n>l. (2.56) 

Vise versa, if we have a basis of the form (2.56), the corresponding point of Gr leads to a 
solution of the i^-reduced KP hierarchy. In the next section 2.6 we show that the GKM 
(1.2) with arbitrary polynomial potential V{x) corresponds to the V'(/x)-reduction of the 
KP hierarchy. 

Note that i^'-reduction is equivalent to existance of some admissible basis with the 
property (2.56). However, other admissible basises (including canonical) in this case 
satisfy a weaker condition: 

n 

4>K+n{f^) = f^^ ^i^i 
1=1 

with some constant a^'s. 
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Fermionic formalism. Grassmannian approach is in fact equivalent to the fermionic 
language used in the section 2.3. Let us briefly describe the relation between them. 

In the fermionic interpretation H becomes the one-particle Fock space spanned by ipn- 
The fermionic operator G of the form (2.26) makes linear transformations in H according 
to 

V^n ^ G^PnG-' e H. 

The positive (negative) modes of the current J(/x) (2.26) are the generators of r_,_(r_). 
The matrix Qmn in (2.26) can be expressed through the canonical basis 4'^k\l^) of the 
corresponding W as follows: 

Smn + G^n = / '^iTv^Gwili, 1^), (2.57) 

CO 

gwM = Y.l^''^k\^)- (2-58) 

k=l 

A more invariant expression is 

1 z/~*^ — uT^ 
Gwil^,'^) = rw({ : }). 

/jL — U K 

Eqs.(2.57) and (2.58) are very important for any interpretation in terms of (infinite-genus) 
Riemann surfaces, in particular, for identification of GKM with appropriate Liouville-like 
models of 2d gravity. Explicit formulae for Qy^ can be easily derived, at least, in the case 
of V = Vk- 

T -functions of reduced hierarchies. Let us note that in the fermionic language it is 
trivial to understand that the X-reduction condition leads to the r-function r^^^ almost 
independent of times T^k, i-^- the property (1-21). Indeed, (2.56)-(2.58) imply that the 
condition of such reduction can be transformed into the matrix {Qmn} '■ 

where E is the shift matrix, {Emn) = 5m+i,n- On the same "infinite-matrix language" the 
elements of F-i- can be parameterized as 7+ = exp(J2TkE''), 7_ = ex]){J2TkE~'^), i.e. the 
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currents Jk from eq.(2.26) act as . Now it is evident that exponent of currents Jk with 
A; = mod K can be pulled to the right vacuum and canceled, eliminating the dependence 
of the corresponding times (up to the exponential factor like that in eq.(2.54) — this is 
due to the freedom to insert the exponential 7_ = exp(X] T/c-^'"^) between G and 7+). 

Now let us discuss the freedom in the definition of the r-function from the viewpoint 
of Hirota equation. As it was pointed out in (2.54) there is a freedom to multiply the 
r-function by a linear exponential of times. Another possible modification is dealing with 
other components Gr„ of Grassmannian. It corresponds to the shift of the first basis 
vector or, equivalently, multiplying it by /i". All these can be trivially understood in 
terms of Hirota equation in Miwa parameterization (2.42). Indeed, it is invariant with 
respect to multiplying by the product Hi /(/^i); where /(/i) is arbitrary Laurent series. If 
the first term of this series ~ one can expand /(/ij) — > Q^viY^k ^M^Tk} and reproduce 
correct factor in eq.(2.54). If the first term behaves as /i", one should deal with Gr^. Thus, 
in Miwa variables we readily describe all full freedom in the definition of the r-function. 

GKM and specific points in Grassmannian. To conclude this subsection let us note 
that the problem of correspondence between matrix models and the points of Grassman- 



nian has been already addressed in ref.p3[. That paper examined the implications of 



£_i -constraint in the double-scaling limit of if-matrix model and concluded that the 
corresponding point of Grassmannian is defined by the set of K linear independent vec- 

XK+l 

tors {A^^i}, where operator A for V = essentially coincides with ours (see (2.19), 

K -\- \ 

(2.20)). Moreover, from our study of GKM we obtain straightforwardly the explicit form 
of all basis vectors for any potentials: (pi = $i as given in eq.(2.18). Note also that in |2^ 
a shift of time variables analogous to ours in (1.17) was introduced. 

As to ref. |2^, we would explain the appearance of the Airy functions (in the case of 
K = 2 model) from constraint (1-20) in the following way. If we put all = except 
for Ti this constraint can be satisfied only if the r-function is zero. This means that the 
relevant point of Grassmannian is singular. Instead, we can put T3 = c 7^ keeping all 
other Tk with k > 5 equal to zero. In this case we have 

3 dloq T 1 9 
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Choosing c = —2/3, so that T3 = T3 + 2/3 = 0, and taking the derivative with respect to 
Ti and using the relation 



<9^1og T 

<T,) = , 

where u{Ti) is the KdV potential, we obtain 

2«(ri) = Ti . 

According to the ideology of the inverse scattering method in the case of 2-reduction we 
should consider the Schrodinger equation for this potential: 



I.e. 



Solving this equation, we obtain Airy function 



+ Ti* = ^i^m . (2.59) 



^(/ilT^i) It,>3=o = '^MAi{~T,+ii'') = vi/o(/i) I e-^'/3+(M^-T,)x^^ _ 

The wave function \E' is nothing but the Baker- Akhiezer function. \E'o is adjusted from 
normalization requirment ^(/i)e~S^"'^" — 1 as — oo. In general \Ef can be expressed 
through r-function as Q 

vI/(/.|T„,) = e^^-z^-li^^i^ (2.60) 

(see sect. 3. 3 below for more details). In eq.(2.59) we have \Ef-function of the form (2.60) 
with Tfc = 0, A; > 3. The admissible basis for the corresponding point of Gr^ can be 
constructed by means of the \l/-function as follows pO||: 



0fc(/^) 



all Tfe=0 



The basis vectors, constructed according to this expression, obviously coincide with those 
which were obtained in sect. 2. 2 (eq.(2.6)) for the particular case of V(X) = ^X^. 
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For generic K the role of (2.59) is played by a K-i\i order differential equations [L^^! ~ 
11^'^) and the result is the Airy function of level K as defined in (2.6). 

2.6 GKM and reductions of KP-hierarchy 

We prove here that whenever potential V{X) in GKM is a homogeneous polynomial of 
degree X + 1 , i.e. V{X) = canst ■ X^+^), the partition function Z{^>[M], which is a KP 
T-function, in fact can be treated as a r-function r^^^ of X-reduced KP-hierarchy. 

For this purpose let us return to the section 2.2 and note that because of (2.7) -Fi(A) 
satisfies the Ward identity (1.23): 

[V{d/dX) - A]Fi(A) = 0. (2.61) 

If the potential V{X) is a polynomial of finite degree K + 1 , then (2.6) may be used to 
express Fk-\-i in the form of linear combination of Fj's with i < K. Namely, if V{X) — 
-ES'V,X^ Vk+1 = eq.(2.6) implies, that 

K 

Fk+1 = - E i^i^i + ^^1- (2-62) 

1=1 

Clearly, all the terms on the r.h.s. of (2.62) except for the last one AFi will drop out 
of the determinant (2.9), i.e. in this case F^+i may be defined as AFi rather than by 
eq.(2.62). Obviously in the same way any Fx+m may be substituted by XFm , while any 
FnK+m — by X^Frn- In other words, J?-7v[A] is given by eq.(2.9) with the first K functions 
Fi . . . Fk defined by (2.6), while other elements of the basis are given by the recurrent 
relation: 



FnK+^ ~ A"F^. (2.63) 

Note that this is true when V{X) is any potential of degree K +1, not obligatory homo- 
geneous. However, we should recall, that A = V'(//) and /i rather than A is the proper 
parameter to deal with in Grassmannian picture. Therefore eq.(2.63) implies the appear- 
ance of KP-reduction, associated with the function V'(/i). If V is further restricted to be 
XK+i 

V{X) — — , A = jj,^, and (2.63) acquires the form: 
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FnK+m ~ A"F^ = (2.64) 

which is characteristic of conventional i^-reduction. 

As we aheady explained, i^'-reduction implies the relation (1.21) with some factor of 
the form (1.22). However, even this factor is absent in the case of GKM (1.10). For 
particular case of V{X) = const ■ {K = 2) this was exhaustively proved in [T^, using 
the properties of symplectic structure on (the model of) the universal module space. Since 
interpretation of entire GKM in such terms is not yet available, the analog of such proof 
for generic K is still lacking. 

2.7 On r„/^-independence of Z^^^ 

Reformulation of the problem. First of all, let us choose our variables {/ii} in such a way 
that the only non-vanishing times Tn = l^i^ are those with n = mod K. To do 
this it is enough to have only K finite parameters /ii, . . . , ^k-, which are essentially K-i\i 
order roots of unity: 

^- = ^e^^ e = e^"/^, j = 1...K. (2.65) 
For such choice of variables the formula (1.21), 



T {Ti . . . Tk-1, Tk, Tk+1 .... T2K-1, T2K, T2K+1 ■ ■ ■) — 
= ei:„««i^^"KrW(Ti . . . Tk-i, 0, T^+i .... T2K-1, 0, T^k+i ■ ■ •), (2-66) 

turns into: 



rW(o ...0,Tk,0...0,T2k,0...) = eS„'^"^^"ifr{^}(0) = e^n ^'^""''r^^'}(0). (2.67) 

This relation is valid for any r-function r^^^ of i^-reduced KP hierarchy. Our purpose 
is to prove that for the particular example of such t^^^ , namely, for partition function 
Z^^> of GKM, associated with potential of the form V(X) = X^+^/{K + 1) , all the 



= d log Z^^^dT^K = 0, 
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(2.68) 



or, equivalently, the function in the r.h.s. of (2.67) is independent of /x. In other words, 
we need to prove that 



e{^>(/.)^zW[/.,]|p.e5) = l. (2.69) 
The l.h.s. of this relation can be evaluated with the help of eq.(2.16), 

{K} ^ dct(,,)'l'^(/ij) 

A(Af) ■ 

We emphasize that in order to prove TnK- independence of Z^^^ it is enough to examine 
the determinant oi K x K matrix in (2.16) with entries 

Hp^) = (^(^)|;^(^)"') ^ ^iK}^i = /^^-'e«(/.), (2.70) 

e«(/^) = EeiV"^^^^U« = l (2.71) 

Q>0 

(of course, all the quantities s, A (see (2.15), (2.18)-(2.20)) and e*^*^ depend on the form 
of potential V and thus on K) . 

Straightforward calculation {K — 2). We begin our proof that ^^^^{fJ>) — 1 from direct 
calculation of this quantity. If we substitute the expansion (2.71) with fj, 's defined in 
(2.65) into (2.16) the result reads: 



K 



1 detr,-,-^£(-°^'+^-^)^' 



The remaining determinant at the r.h.s. is actually vanishing unless J^iCa — ^ mod K, 
so that (2.68) is in fact equivalent to the set of identities 

for coefficients of series expansions of modified X-level Airy functions 

$f >(;,) = Vg-^^ J e- W+'^^V-irfx. (2.74) 
For example, for X = 2 identities (2.73) look like 
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a+/3=2n 

Since for X = 2 



E = (2.75) 

a+/3=2n 



1 r(3A; + l/2) 1 



E 



$2(/x) 



9^'(2A;)! r(l/2) /x^^ ' 
■^e~^^^ J dx- X 6"^+^''^ = (i2;(2; + /x) exp j-^;^ - ^ 



3/^3/2 



^ 1 6A; + ir(3A; + l/2) 1 



k 



9'=(2A;)!6A;- 1 r(l/2) //^fc ' 



i.e 



,(1) - 1 r(3a + 1/2) (2)^_6a + l 

9"(2q;)! r(l/2) ' " Ga-l " ■ 

A more explicit form of (2.75) is 

^ ^WQ^Il ^^^" - V2)r(3/3 - 1/2) = -4[r(l/2)]^^.,o , (2.76) 

a,l3>0 

which is indeed a vahd F-function identity. 

This calculation, though absolutely straightforward, is hard to repeat for generic K. 
Still it is useful for illustrative purposes and we included it into this section. 

The proof. A much easier approach is to prove that 

/^|;e^^>(/^) = 0. (2.77) 
Since it is obvious from (2.72) that at least C,'^^^ — 1 as /i — > oo, this would provide a 
complete proof of (2.69). In order to prove (2.78) it is enough to act with iid/d/j, upon 



e{^>(//)~det(,,)£^^e«(/x£^) (2.78) 

and then make use of the fact, that fi-^e^^\iJ,e^) can be decomposed into linear combina- 
tion of e^^\fj,e^) . . .e^^\fj,e^). This decomposition, which is, of course, the implication of 
reduction, can be worked out from the relation (2.20), 
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Since 



(2.79) 



d K-l 

(see also lESl), we have: 



(/x|-)e«(/^^') = i ^^l - ire>^+')e«(/x£^) + ir£V^+ie(^+i)(^£^). (2.81) 
This should be supplemented by the condition 



e(^+i) = e(i) (2.82) 
(equation Ae*^^^ = ji^e^^^ = A^j^^e^^^ = s(A) f — j s{X) ^e^^^ is nothing but the equation 
for the level-K Airy function F^^^ = s^^e^^^ = / e~^^K+r+i''^''dx). 



Let us begin our study of (2.78) from the familiar example of = 2. First of all, one 
can check up (2.82) with the help of explicit expressions (2.76) for eW. Then 

- (/i|^)(e«(/.)e(^)(-/.) + e(^)(/.)e«(-/.)) = 
= [(i - 2/.^)e«(/.) + 2/.V^)(/.)]e(^)(-/.) + e«(/.)[-2/.^e«(-/.) + (-1 + 2^,')e^^\-^^)] + 

+ [2/iVi)(/i) + (-^ - 2/.3)e(2)(^)]e«(-^) + e(^)(^)[(i + 2^^')e^'\-^^) - 2f^'e^'\-f^)], 

and the r.h.s. is identically zero. 

If we return to the case of arbitrary K , note that after the substitution of (2.79) into 
the l.h.s. of (2.78) we obtain a sum of K determinants like (2.79) with e^'^e^''''\fie^) in the 
row i' substituted by (2.82): 



Next, note that the last item in the r.h.s. coincides with the entry e'^'^' ^^^^ e'^'^' i^iie^) of 
the next line of the same determinant, and thus can be eliminated. Moreover, the first 
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K -\-l — 2i' 

item at the r.h.s. implies just that the i'-th Une is multiphed by a factor of ^ 

K — 2i' 

and the effect of this cancels in the sum over i' : J2i/=i 7: = 0. Therefore, we 

conclude that 

f'^^^'^Hf^) ~ /x|^detfe)£^^e»(/.£^) ^ t,-^V = J2 V, , (2.83) 

where Pj/ is just the same determinant T) , only in the i'-th line di'j = e^'^e^'^'\iJie^) is 
substituted by e^diij = e'^^'^^^^ e^'^'\jj,£^). The sum of such (for any original V) is 
identical zero as a consequence of the identity J2f=i — 0. 

Examples: 
K^2: 



du di2 

= det I I = dnd22 - c?i2C?2i; 

d21 d22 

-dn i-Ydu \ 

(^21 d22 



Vi = det 



"<^ll<^22 ~ C?12<^21) 



(ill di2 

= det = dud22 + c?i2C?2i, 



-d21 {-fd; 



22 



and 



Pi + P2 = 0. 



X = 3 : £ = e2'^^/^ 



V = det 



dii 


di2 


di3 


d2i 


d22 


d23 


dsi 


dz2 


d33 



dlld22d33 + 0?120^23C^31 + 0^13^^320^21 ~ C^12C^2lC^33 ~ C^13C^3lC^22 ~ 0^23^^320^11) 
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Vi = det 



C^21 (^22 C^23 
C^31 C^32 C^33 



SC^ll'^22'^33 + £^'^12'^23'^31 + £^^C?13C?32C?21 ~ £^(^12(^21(^33 ~ £'^C^13'^3lC?22 ~ £'i23C?32'^ll, 



X>o = det 



V 



dn 


di2 


di3 


ed2i 


e'^d22 


£^d23 


dsi 


d32 


d33 J 



£^(^11^22(^33 + £^(^12(^23(^31 + £(^13(^32(^21 ~ £(^12(^21(^33 ~ £^(^13(^31(^22 ~ £^(^23(^32(^11, 



^ (ill di2 dis ^ 



©3 = det 



(^21 (^22 (^23 
-^'^33 



\^ edsi s2d32 £ (^3 

= £^(^11^22(^33 + £(^12(^23(^31 + £^(^13(^32(^21 ~ £^(^12(^21(^33 ~ £(^13(^31(^22 ~ £^(^23(^32(^11, 



and 



Pi + P2 + 2^3 = 0. 

These two examples are enough to illustrate the general phenomenon. The thing is 
that determinant V is an algebraic sum of products di^j^ ■ ■ ■ di^^jj^. Let us pick up one of 
these items. It appears in T>i> with the coefficient e^' , where j' is defined as the second 
subscript of the element di'j/ in this product. The map i' — > / depends on the particular 
item we choose, but when we sum over all i' this is the same as to sum over all f , since 
every subscript appears in the product once and only once. Therefore in the sum J2i' T^i' 
every item appears with the universal coefficient X^^-z e^' — , i.e. T^i' = 0. 

This completes our proof of (2.78) and thus of (2.68). 
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Remark. To avoid possible confusion, let us note that it is important that Airy- 
functions are analytic functions of /i , but not of A = iJ,^{e.g., e^^\fx) are expanded in 
series with integer powers of jj,^^^^^^ = \^^^^^^). This is what makes the proof slightly 
non-trivial. Otherwise, if one suggested that the entries of the matrix in (2.79) are analytic 
functions of A , the A-derivative of this determinant would be simply a combination of 
determinants with coincident rows and thus vanishing. However, non-analiticity in A 
requires one to be more accurate: X^^^ takes different values fxe^ when appearing in 
different places, and a detailed calculation, as given above, is necessary. 

3 Universal £_i-constraint and string equation 
3.1 Motivations 

Since it was proved in sect.2 that Z^^^M] is a r-function t^^^ of the i^-reduced KP 
hierarchy, in order to define it completely {i.e. to fix the point of Grassmannian) it is 
enough to deduce somehow a single additional constraint 

jr^zW = 0, (3.1) 

with 

CiV = W^^l = ^ E in + K)Tn+K^/^Tn-^/^T^ + ^ ^ aTabU. (3.2) 

n>l a~^b=K 
nj^O mod K a,6> 1 

According to arguments of ref.[|l| such constraint (3.1), when imposed on Z^^^ = t^^\ 
implies the entire set of PVi^-algebra constraints in the form of (1.29). Another reason for 
the study of C-i- constraint is that it is much simpler than any other ones: C-i is the 
only operator of interest which does not contain double- and higher- order T-derivatives 
(to be exact, there is one more such generator: Cq). Among other things, this means that 
it is easier to formulate the universal {i.e. for any potential V{X)) £_i-constraint, than 
any other corollary of (universal!) Ward-identity (1.23). 

This section contains formulation and derivation of iZl^^-constraint for GKM with 
arbitrary potential V{X) , making use of explicit formulas, derived in sect.2. Namely, we 
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are going to exploit the fact that the functions in (2.10) are not independent, but 

are rather related by the action of d/dX- operator: see (2.20). The proof itself is described 
in sect. 3. 2 and sect. 3. 3 contains additional remarks about string equation. We emphasize, 
that not only the C^i- constraint is valid for any V{X) (with the only restriction that V(x) 
grows faster than a; as a; ^ oo), it is just the same for any size of the matrices. Just 
like the property of integrability (i.e. that Z^"^^ = r^"^^), this constraint is not sensitive 
to , and this makes the entire construction behaving smoothly in continuum limit, as 
N oo. 

3.2 Direct evaluation of jC-iZ 



It is well known pi] , |28| , that £_i-constraint is closely related to the action of operator 



Therefore it is natural to examine, how this operator acts on 



Zi^^m-'-^Ilsi,.), (3.4) 



s(^) = (V"(/i))^/'e^(^)-'^^'(^), (3.5) 



l>,(/i) = F,{X) = id/dXy-'F,iX), X = V'ifi). 
First of all, if Z^^^ is considered as a function of T-variables, 

In the particular case of V(X) = Vk{X) = — , eq.(3.6) turns into 

i\ + 1 

1 ^ J_^{K} ^ V T ^ ^^^""^ = 

Zi^} ''dAtr KZi^} ^ ''M"+^ dTn 

71^0 mod K 

' E aT.6T, + AUW (3.7) 



^^^'^ 1 2 dT, 

a,b>l 
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Wc made use of the definition (3.2) of -operator and tiie fact that Z^^'^^ is independent 
of all (since it is a r^^^ -function). Time- variables are defined by Miwa transformation 
Tn — ^Tr M~^. On the other hand, if we apply (3.3) to explicit formula (3.4), we obtain: 



1 -Tr^Z{^> 



^-^ + 22. V"(/.OV"(/.,) Pi^ - Pij + ''dA,r 

or, in the case of V{X) = X^+^/{K + 1) , 



^rr^Z^^> = -TrM +\ Y.^ <^T,bT, + Tr-^^ log det F,(A,). (3.9) 

a,6>l 

Combination of (3.7) and (3.9) implies, that 



The r.h.s. here is practically independent of K , and this may be used, together with 
eqs.(3.6) and (3.8) in order to suggest the formula for the universal operator jC^^^ : the 
idea is to preserve the relation (3.10): 

= - log Zj^^ + TrM - Tr-^^ log det F,(A,). (3.11) 

Here 



{V} ^ ^ 1 d 

~^ ^V"(M)M"+i^(9T„ 

and this turns into (3.2) when V(X) = X^+^/ {K + 1) (note that the items with i — j are 
included into the sum at the r.h.s. in (3.12)). 

So, in order to prove the -constraint, one should prove that the r.h.s. of (3.9) 
vanishes, 
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log = TrM - Tr-^^ log det F,(A,), (3.13) 
for — rj^^, defined by (3.4). The l.h.s. may be represented as residue of the ratio 

-.^^%t£ZM = Jr.og.rW. (3.14) 

However, if expressed through Miwa coordinates, the r-function in the numerator is given 
by the same formula (3.4) with one extra parameter /x , i.e. is in fact equal to rj^^i . 
This idea is almost enough to deduce (3.13). Since (3.13) is valid for any value of N, it is 
reasonable to begin with an illustrative example of N — 1. Then (A = V'(//)) 

rP^(7; + //-7^)=rf = 



[F{Xi)dF{X)/dX - F{X)dF{Xi)/dXi] 
/i — /ii 



(1- (ll 

The function 



rPil^i] ■ [-dlogF{Xi)/dXi + dlogF{X) / dX]. 



(3.15) 



r V"" 
F{X) = dx e-^(^^+^^ - e^^^^-^^'^^\V'\ii)]-^'\l + ^(-J— )}. (3.16) 

If V(/i) grows with as ^ oo , then V"" /(y"Y ~ //~" , and for our purposes it is 
enough to have n > 1 , so that in the braces at the r.h.s. stands {1 + o{l/ ^)}{^ ■ o{fi) — > 
as fi oo). Then numerator at the r.h.s. of (3.15) is ~ 1 + o(l//i), while the second 
item in square brackets behaves as dlogF[X)/dX ~ /x(l + o(l//x)). Combining all this 
with eq.(3.14), we obtain: 

-^logr^ = res, ( ^ + "^^^^^ [-dlogF{X,)/dX, + /x(l + o(l//x))] 1 ^ 

= Hi-dlogF{Xi)/dXi. (3.17) 

Thus (3.13) is proved for the particular case of A?" = 1. 
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The proof is literally the same for any , only instead of a relatively simple expression 
in square brackets at the r.h.s. of (3.15) one has the ratio: 



det 



F(Ai) 9F(Ai)/9Ai 

F(A2) dF{X2)/dX2 

F(Aiv) 9F(Aiv)/9Aiv 

F{\) dF{X)/dX 



a^-iF(Ai)/9Af-^ a^F(Ai)/aA^ 



d^-^F{X2)/dX^-^ d^F{X2)/dX 

d''-'F{Xr,)/dX^-' 9^F(Aiv)/9Aj^ 
9^-iF(A)/9A^-i a^F(A)/9A^ 



(3.18) 



over 



det 



F(Ai) dF{X,)/dX, 
F{X2) dF{X2)/dX2 



d^-^F{Xi)/dX^-^ 
d^-'F{X2)/dX^-^ 



F{Xn) dF{XN)/dXN 
which is in fact equal to 



d''-^F{XN)/dX 



N 



dF/dX , d^F/dX^ 
We used here the estimates — - — = //(I + o(l///j) , 



(3.19) 



F(A)/.^ |[1 + 0(1//^)] - -^[Tr^^Xog det F,(A,)] • [1 + 0(l//x)]| . (3.20) 



o(l//x))]=/x2(l + o(l///)), 



9^F/9A 



N 



dX' 



fjL (l + o(l//x)), which allow us to pick up 



only the contributions with d^F{X)/dX^ and d^-^F{X)/dX^-^ to (3.18) — all other are 
lower order in as — > oo. The A^-th derivative in (3.18) is multiplied by determinant 
of N X N matrix, which is exactly (3.19), while the {N — l)-th derivative — by 



det 



F{X^) ^F{X^)/^X^ 
F{X2) dF{X2)/dX2 



9^-2F(Ai)/9Af-' 9^F(Ai)/9Af 
9^-2F(A2)/aA^-' d'^F{X2)/dX^ 



F{Xn) dF{XN)/dXN ... d''-'F{XN)/dX^-' d''F{XN)/dX^ 



(3.21) 



d 

This is, however, nothing but Tr— — of the logarithm of (3.19), and for this to be true 
it is absolutely essential, that Fj ~ {d/dXy~^Fi. 
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With the estimate (3.20) we obtain from (3.14): 



l + o(l//x) 



N Q 

J2 N - Tr—- log det F,{Xj). (3.22) 

j=l '^^Hr 



This completes the proof of eq.(3.13) and thus of the universal vCl^-'^-constraint: 

^m^m _ Q_ ^323) 

We proved the crucial eq. (3.13) by the direct calculation. Now we would hke to describe 
an alternative proof which originates from representation of r-function through the current 
correlators described in sect. 2. 3. This approach may be useful for evaluation of the higher 
T-derivatives (see sect. 3. 4 below). Namely, from eqs.(2.35) and (2.38) it is obvious that 
T^[t\M] = r{t, {Tn}) depends on t and as follows: r{t, {r„}) = r{Ti-Nt, ^2-|A^^^ -) 
i.e. in eq.(2.38) 

f"' 8 

T^[t\M] = exp(-7V^ )r(r) = Y.t^Vn{d)r{T) (3.24) 

where Vn{d) are the Schur polynomials and symbol d represent the vector with compo- 
nents {-Nd/dTi, -\NdldT2 , ...). Therefore from eqs.(2.38), (2.39) and (3.24) one can 
deduce 

-N-^j{T)^d,TS[t\M] |^^^ = -Ar^;,,.r(r) + ^lima,det4(/.,-,i) , 
where (f)i{iij,t) are defined by eq.(2.39). From this expression it follows that 

d 



T{T) = Y.l^r-r{T)-f{T) (3.25) 
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where f{T) denotes some new r-function which is obtained from the given t(T) by the 
shift of the last row in the determinant: 0Ar(/ij) 0Ar+i(/ij). Eq.(3.25) is nothing but 
eq.(3.13)g. 

3.3 Universal string equation 

The string equation, as imphed by (3.23), is by definition: 

If we substitute exphcit expression (3.12) of the £_i-operator the string equation acquires 
the form: 

^ 1 d'^logr 

nl^i V"(M)M"+iJ9T\97; " ■ ^ ' 

Here we introduced a conventional notation u = d^logr / dTf and also defined new Tq- 
and T_i-derivatives, so that dlogr/dTo = , dlogr / dT_i = T^. The derivation of (3.27) 
involves taking Ti-derivatives of the objects 

In order to get some impression about OTn/dTi, let us imagine, that V(M) is some 
polynomial of degree Q + 1 , so that [V"(M)]-i ~ M^-'^{1 + viQ-^ + ^2^"^ + ...). Then 
% ~ {Q + n)TQ+n + vi{Q + n + l)TQ+n+i + V2{Q + n + 2)TQ+n+2 + ■■■ ■ Therefore whenever 
Q + n > 1 {i.e. Q > 0, as it is already necessary for the derivation of (3.23)) 

dTn/dTi = 0, for n > 1. (3.29) 
In order to derive (3.27) one also needs Ti- derivative of the second item in (3.12), 



■^Let us emphasize that exphcit form of r.h.s. of (3.25) does not depend on the choice of basis. Unhke 
this, the manifest form of derivatives over T™, with m > 1, should deal with concrete choice of basis. It 
is noteworthy to remark that, for m < K, this manifest form coincides for the canonical basis and our 
basis (2.18). 
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Under the same assumptions about V(/i) the second ratio in this sum is a polynomial in 
/i's of degree Q — 2 . The only contribution to (3.30), which contains pure factors or 
, i ^ — ^ (1^ + l^U T.r_„ othe.= a.e expanded a= bUinea. 

series in T^T^ with m, n > 1. Therefore, Ti-derivative of (3.30) (if Q > 2) is just Ti and 
is described by the n = term in (3.27). 

So, we derived the universal string equation (3.27) in the form: 

n>— 1 ^ " 

If potential V(X) = const-X-^+i , the variables T„ = , and we return to familiar 

relation p4| : 



— 2^ n + i^ T^+i^ = M . 3.32 

n^O inod K 

Note that even if we restore the factor (1.22), describing the maximal possible Tnx- 
dependence of t^^^ it drops out of the string equation (3.32) after taking the d/dTi- 
derivative. 

Our next task is reformulation of (3.31) in terms of pseudo-differential operators and 
in the form of bilinear relation. If 



L = d + Y.u.^.d'^ id^QjT^Q^)^ (3-33) 



then [33 



d^logr 



in-i . (3.34) 



OT.dT^ 

For n = —1, 0, 1 we have (L^^)_i = 1 , {L^)-i = 0, {L)_i = U2 = u in accordance with 
(3.27). K-reduction is associated with the condition 



(L^)_ = (3.35) 

(this guarantees, that items with n = mod K do not appear in (3.32)). Rewritten in 
these terms, eq.(3.31) turns into: 
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rry^Tp^ E^M-"-(L")_,=«. (3.36) 



Recall now, that in the framework of the dressing formalism |33| L = KdK ^, L" 



Kd^K~^. The eigenfunctions of L can be defined as {Ke^^), where brackets denote that 
operator K acts on e^^ {i.e. [Ke^^) is a function of z , not an operator). Baker- Akhiezer 
function in these terms is given by |Q \1/(2;|T„) = e^^"^" (Ke^^), while its conjugate 
^*(2|T„) = e-S^"^"[(/s:-i)*e-^^] . These definitions are useful for us, since 



(L")_i = res,z"^*(2)^(2). (3.37) 

We assume that contour integral over z , implicit in the definition of res^, is around zero. 
If (3.37) is substituted into (3.36), we obtain: 

M ^f*{z)^{z) 

u = —reSzTr————— — — . 3.38 

V"{M)z{zI-M) ^ ' 

Since ^{z\Tn) = eE li^^,i-!M and ^*(z|TO = e-E^-^-li^^^I^M, the 
product '$*{z)'${z) ~ 1 + 0{l/z) (as z oo) and the contour integral over z , when 
pulled to infinity, picks up the contributions from eigenvalues of M . Finally we obtain our 
universal string equation in the form of a bilinear relation for Baker- Akhiezer functions: 

^ ^*(Ui)^(Ui) 

Comparison of eqs.(2.46) and (3.39) enables us to derive a useful representation for 
M-function: 

Note that eqs.(2.45) and (2.46) which follow from the general formula (2.44), provide the 
connection between residues of the Baker-Akhiezer function and time derivatives of the 
r-function without reference to the pseudo-differential operator language. Indeed, it can 
be easily understood if one substitutes eq.(3.25) into eq.(2.45). 
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3.4 Discussion 

In this section 3 we derived the universal £l^-^-constraint (1.25) on partition function 
Z^'^^ Its simplest implication, the universal string equation, was considered in sect. 3. 3. 
String equation is nothing but Ti-derivative of this £_i-constraint. Remarkably enough, 
it appears to be much simpler than the constraint itself (compare (3.12) and (3.40)). 

However, string equation is not the only implication of (1.23). Another one should be 
the entire tower (1.29) of Virasoro and W-constraints imposed on Z'^^J'. In this subsection 
we discuss three different ways to derive these remaining relations for any V(X) in the 
GKM. Being technically different, they emphasize different (though related) properties of 
partition function, and therefore each of the three deserves careful investigation, which 
is, however, beyond the scope of this paper. 

KP-hierarchy approach. This one is the most popular in the literature (see ||l], Q). 
The idea is that £_i-constraint, when imposed on KP r-function, automatically implies 
all other constraints. In other words, £_i-constraint is enough to fix the point in Grass- 
mannian with which the KP r-function is associated and, thus, it fixes this r-function 
completely. It would be interesting to find out explicit dependence of the point in Grass- 
mannian on the potential V{X) in GKM. Even more interesting would be any alternative 
description of the entire subset Uoo in Grassmannian, associated with GKM with arbi- 
trary V{X), and any reasonable parametrization of Uoo, which would be an alternative 
parametrization of the space of potentials. Since V(X) can be changed smoothly, Uoo 
should be a kind of a manifold, lying at infinity of the universal module space (if the 
latter is embedded into Sato's Grassmannian), and surely possess some amusing proper- 
ties. Unfortunately, the already existing papers on the subject (of course, they concern 
conventional (multi) matrix models rather than GKM) either emphasize implications of 
i^T-reduction, like [Q, or rely upon the formalism of pseudo-differential operators, like 
im, and thus need be translated into Grassmannian language. 

A separate question in the framework of this approach is why should the constraints 
be always associated with Woo-algebra (as they actually are). Technically this is more 
or less obvious whenever pseudo-differential operators are used. In KP- Grassmannian 
language, the crucial ingredient should be the Woo-covariance of KP-hierarchy, discovered 
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in [pSj p6| , |37| . The constraints of interest involve operators from some Borel subalgebra 
of this Woo, and the choice of subalgebra depends on the choice of potential. 

Reductions of GKM to (double-scaling limit of) multimatrix models, i.e. the what 
happens if V{X) ~ X^~^^, in Grassmannian language should be attributed to intersections 
of Uoo and submanifolds Gr^^^ in Sato's Grassmannian, associated with /^-reductions. 
Somehow at these points a closed Zamolodchikov's WA-subalgebra emerges from entire 
Woo- It is certainly interesting to see, how this happens and how the Lie algebra structure 
is broken. It is also unclear, whether such phenomenon takes place only at Uoo H Gr^^^ 
or everywhere in Gr^^\ 

One more important question is what is the relation between the just discussed Woo- 
algebra of [^, ^] and another Woo, which is presumably relevant in c = 1 models, and 
is naively generated by operators, which involve finite-differences instead of derivatives in 
the free-field representation. We use this chance to note, that the problem of how c = 1 
models are included into framework of GKM is still open, and the invariant description 
of the Uoo subset in Grassmannian would be also helpful for its resolution. 

Straightforward derivation of C- and W- constraints. This approach is just the straight- 
forward generalization of the what was done in sect. 3. 2 in the derivation of £l^^-constraint 
from the knowledge of explicit expression (2.2) for Z'^^^. We sketch here several steps of 
the derivation for the Virasoro case. The idea is to apply the operator TrA™"^^g^ to 
(2.2). On one hand this is equal to 

„ji VV"(A^)A^"+V aT„ • ' ' 

on the other hand, it can be explicitly evaluated. Expression (3.41) is very close to 

£i^>Zii^}. The main difference is that the sum in (3.41) goes from n = 1 and for small 
/ (V'(M))'"+i \ 

enough n Tr \ yr/^j\^-^j\^n+i j contains powers of positive powers of M, which can not be 
expressed through T-variables. For V ~ X^~^^ this happens for n < Km. In order to get 
rid of these positive powers of M one should use the analogues of eq. (3.14), saying that 

d 



log Z^^^ ^Tr + ... . (3.42) 
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This is the origin of conventional d'^ / dTadTb-terms in operators C^^\ As to generahzation 
of the shift (1.17), it results in additional contribution of the form 



in the definition of C^^^ - operators. All extra terms cancel, giving rise to the proper 
universal constraints of the form Cj^^Z^"^^ = 0, m > —1. 

W-constraints can be (at least in principle) derived in the same manner, though actual 
calculations are increasingly sophisticated. 

Implications of Ward identities. This third approach exploits the fundamental Ward- 
identity (1.23). Technically its main difference from the previous approach is that the 
operators TrV'(^^), non-linear in d/dAtj., are involved. Conceptually this approach is 
different, since it does not exploit explicitly any integrable features of partition function. 
We shall discuss this method in more details in the next section 4. 



4 From Ward identities to VF- constraints 
4.1 General discussion 

This section is devoted to the derivation (unfortunately, incomplete) of the entire set 
of Virasoro and W- constraints in GKM. The role of such constraints in the study of 
any matrix model is two-fold. First of all, they can be considered as complete set of 
differential equations, which specify partition function of the model as a function of time- 
variables. Then this set of equations implies, among other things, that the solution is KP 
r-function (and sometimes a reduced r'>^^J'-function). The first application, discussed in 
details in sect. 3 above, is reasonable if one already knows about the integrable structure. 
Then partition function is identified with some (reduced) KP r-function {i.e. evaluated 
at some point of Grassmannian) , and the role of Virasoro constraint is to specify this 
r-function {i.e. fix the point in Grassmannian). In this second type of circumstances it is 
enough to have a single constraint, namely C^\Z = , — all other constraints, if imposed 
on KP r-function or t^^^ , are deducible corollaries of this one [ffl, pOl- With identification 



of with a r-function in sect. 2 and the proof of universal constraint in sect. 3 
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we exhausted this hne of reasoning. The subject of this section is to concentrate instead 
on another approach and ignore almost all what we already studied about integrable 
structure of GKM. Then the constraints can be considered as a complete set of differential 
equations, which specify the partition function of the model as a function of time- variables. 
Among other things, the set of equations implies that the solution is KP r-function (and 
sometimes a reduced r'^-'^^-function). To be more concrete, we shall investigate direct 
corollaries of Ward identity (1.23), 

{Tr e(A)[V '{d/dAtr) - A]}^^[A] = 0. (4.1) 

We shall, however, discuss only specific potentials, V{X) = const ■ X^~^^ , and use explic- 
itly the fact that partition function is independent of all TnK- 

The main problem with the implications of Ward identity (4.1) is that they acquire the 
form of conventional Virasoro or W-constraints only in the limit of = oo. The reliable 
results from our point of view, should, however, be A^-independent. But as soon as A^ 
is finite, the complete set of independent Ward identities (4.1) is also finite. Remarkably 
enough they can still be expressed through generators of W-algebras, but the constraints 
arise as a finite number of vanishing conditions for infinite linear combinations of W- 
operators, acting on partition function. As A^ — oo , the number of such conditions tends 
to infinity, implying that every item in linear combinations vanishes by itself. These items 
have exactly the form of conventional Virasoro constraints in the particular case of K = 2, 
while for K = 3 they rather look like eqs.(1.26), and generalization of (1.26) for A" > 3 
is more or less obvious. In any case the honest statement is that any solution to the W 
-constraints in the conventional form (1.29) does satisfy the equations, which follow from 
(4.1), for any A^ (to make A^ finite one should take all but the first A^ eigenvalues of matrix 
M to infinity, see sect. 2. 2). In this sense transition from finite to infinite A^ is smooth. 
However, inverse can be true, i.e. (4.1) can imply (1.29) only as A^ = oo. Moreover, since 
(1.26), which is actually implied by (4.1) as A^ = oo , is not quite identical to (1.29), one 
needs also to rely upon the (very plausible) assumption that solutions to (1.26) and (1.29) 
are unique (up to inessential constant factors) and thus coincide. 

Such approach has been already applied in [0] to the study of Kontsevich model itself, 
i.e. for the case of A^ = 2 : V2{X) = Our purpose now is to extend this consider- 
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ation to other potentials Vk{X) = X^^^/{K + 1). However, while the structure of the 
answers is very clear (this is obvious from (1.26)), actual calculations are very tedious. 
Therefore we restrict our presentation below only to the first non-trivial case of = 3 
. It is considered in sect. 4. 3. Before, in sect. 4. 2 we reproduce from the derivation of 
Virasoro constraints in the case of i^' = 2. 

4.2 Virasoro constraints in Kontsevich model {K = 2) 

The problem. Our purpose in this section is to prove the identity 

i^K^pT^ - 'p^)^ = ^ E ^nZ tr{e,A-~') (4.2) 

-r OAtr ^ n>-l 

for 



T^'^^{K]= /dX exp(-trXV3 + trAX) = C[VA]exp(-trA3/2)z/f (T^) / m - odd 
J 3 

(4.3) 

with 



C[VA] = det(yA®/ + /(g) VA)"^ (4.4) 

and 



£{2} ^1 kTk—^— + - ^ 

oTam 



k odd a,b>0 ; a,b odd 

Below in this section we omit label {2} and do not indicate explicitly the fact that all 

sums run over only odd times. 

While (4.2) is valid for any size of the matrix A, in the limit of infinitely large A 

(X oo) we can insist that all the quantities 

^To avoid any comparison note that in original version of Q there were wrong signs in the exponential 
in (4.3) and in the shift of times (1-17) (normalization of times are also different ). 
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tr{e,K-^-^) (4.6) 

{e.g. trhP"'^''^ for ep — hP) become algebraically independent, so that eq. (4.2) and (4.1) 
imply that 

CnZ{T} = 0, n > -1 . (4.7) 

Method. Note that ^{A} in (4.3), which we have to differentiate in order to prove 
(4.2), depends only upon eigenvalues {Afc} of the matrix A. Therefore, it is natural to 

consider eq.(4.2) at the diagonal point Ajj = 0, i 7^ j. The only "non-diagonal" piece of 
(4.2) which survives at this point is proportional to 



dAijdAji 



Amn=0, mjtn "^i 



for i ^ j. (4.8) 



Eq.(4.8) is nothing but a familiar formula for the second order correction to Hamiltonian 
eigenvalues in ordinary quantum-mechanical perturbation theory. It can be easily derived 
from the variation of determinant formula: 



Slog{det A) = tvjSA - ^tr{jSAjSA) + ... . (4.9) 
For diagonal Aij — XiSij, but, generically, non-diagonal SAij, this equation gives 



^ SXk _ _ 1 ^ SAjjSAji _ 1 ^ / SAjjSAji ^ 

k j^j XiXj 2.^.\Xi XjJ Xj-Xj 

which proves (4.8). 

Proof. Now we shall turn directly to the proof of (4.2). Since €p is assumed to be a 
function of A, it can be, in fact, treated as a function of eigenvalues A^. After that, (4.2) 
can be rewritten in the following way: 



c{VA)z{t} 



C(VA)ei*'-^'^'z{r} = (4.10) 
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y pfp pfp Z^^pK'^i) ^^ ■ + 



dX dX 



dX V3/ dX 



+ 



+ 



92 



^dAijdAji \3j 



+ 2^6,(A,) (-) 



2\ dtrA^/^dlogC 



dX dXj 



with trA^/^ = Efc Afc/^ and 



1 



First, since 



c=n(\/A,+v^,)-^/2. 



2J 



it is easy to notice that the term with the second derivatives (4.11) can 
written in the desired form: 



^ n>-l i a+fe=2n "-^ 

= i E t.fe(A.)A-»-) E 

n>— 1 a+b=2n " " 

For the first- derivative terms with the help of (4.17), 
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dlogC 



E 



(4.20) 



and 



dKijdKji 



Then for (4.12) we have 



1 ^ dZ 
2^ nVo 



n>0 



^ ri-3/2 , n-3/2 



\/Aj + \/A,- 



j a+b=2(n+l) 



1 



-n-1 



n—1 

E E 6,(A0A-"-^+'^A-"-^/^ - 2 E e,(AOA-"-^ 

i,j a=0 i 



^ E E^p(A.)Ar-^ E (ea7'-^^'-24 



n>— 1 i fe>5„+i,o \ j 



dT2n+k 

dZ dZ 



n>-l 



(4.21) 



(4.22) 



(4.23) 



The remaining part contains terms proportional to Z itself which need a bit more care. 
First, it is easy to notice that two items in (4.13) just cancel each other, so the (4.13) 
gives no contribution to the final result. For (4.14) we have 



Using (4.21) and (4.8) it can be transformed into 



(4.24) 



^Eep(Ai)Ai +2^ep{Xi)— 



Eep(Ai 



1 



(4.25) 
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and this cancels (4.15) (transformed with the help of (4.20), (4.21)). Thus, the only 
contribution is (4.16), which gives 



1 d^C ^ d^Xk dlogC 1 S^C 



E 



t 2(A. - A,) 



\/Xj{\/\j + \/\k) VAi( VAi + v^ifc) . 



(1-5.,)+ 



(4.26) 



+5, 



dHogC ( dlogC V 



Now, using (4.8) and (4.20) we obtain for (4.16): 



+ 



16 



+ iE^-(^^)^^^'(vA,(v^, + A.)^(^. + A)^) 



+ 



- i,j,k ( VAi + VXj) {VXi + VXk) 

_^ 1 ^ £p(At) 

2 i-^j^k -^^ - 



+ 



} 



-3/2 ^ -3/2 



^ E V' + E (^e,(A,) ArU-i + i6,(A,) Ar^/^AT^/^ ) + 



^ iU/"^^'^ (4A,(v^, + VXj)(y/Xi + VXk) 



2{Xi - A,) [ VXjiVXj + VXk) VXjiVXi + VXk) 

^Y.^pi\)K' + E (l^pi\)K'>^j' + ^^p(A.)Ar/^A7^/^) + 



1 



+ 4 E e(A)A-U-/^A-/^ 



(4.27) 



In the last transformation we used the fact that 



4 ^ 



1 1 

VAj ( VAj + VAfe) VAj- {y/Xi + VAjfc) 
1 1 



1 (A. - A,) 



(4.28) 



VAj ( VAj + y/Xk) \fXj{y/Xi + ^^Xk) _ 



+ (J 



fc) 
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Finally, (4.27) can be rewritten as 



16 



n \ i J 



1/2 



> . 



(4.29) 



Now taking together (4.19), (4.23) and (4.29) we obtain our main result: 



c(\/A)z{r} 



C(VA)ei*''^'^'z{r} 



EMe.A-"-^)U E ^-^^^ + 1 E 

n>-l fc>5„+i,o 



(9^ 



a>0,6>0 



+ ^^„,c + iw?-gj^U(r) = o. 



(4.30) 



This completes the derivation of (4.2) and, thus, of Virasoro constraints for the case 
4.3 Example oiK = Z 

In the case of generic K the analogue of the derivation, described in the previous section, 
should involve the following steps. 
— Represent ^[A] as 



^W[A]=^k[A]^^''^(T„), 



(4.31) 



with 



^^^''"^niA.-^e'^^^^^^^^]. (4.32) 



A(A) ' A(A) 

Parameter a is introduced here for the sake of convenience, in fact, a — \. 

— Substitute this ^^^^[A] into (4.1), which in the particular case of Vii-(X) 
looks like 



(¥Ti)' 
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'^^^^^M;/ ~ c.^^]}gK[A]Z^''HTn) = 0. (4.33) 
When d/dh-tr acts on Z{Tn), the following rule is applied: 

|^ = -lg(n.A-)T„,.f. (4.34) 

and the sum at the r.h.s. goes over all positive n 7^ mod K. As to higher-order 

derivatives, 'qJ^ ■• they are defined with the help of relations like (4.8). 

— When all the A-derivatives in (4.1) act on exponent in g[h\ = g\y '(M)] , we get the 

OTriMV '(M) - V(M)) 

term, which is equal to V '( — ) = V (M) = A and cancels A-term 

oAfr 

in (4.1). The next contribution, when all but one of the A-derivatives act on the exponent, 
vanishes. Actually this reflects the fact that there are no W^^^- generators among the final 
W- constraints. 

— Perform a shift of variables 



K 

Tn^Tn^Tn - a.—5n,K+i (4.35) 

77/ 

(this procedure doesn't change d/dTn — >■ d/dTn)- 

— After all these substitutions the l.h.s. of eq.(4.33) acquires the form of an infinite 
series where every item is a product of Tr[€{M)M~P] and a hnear combination of gen- 
erators of W^-algebra, acting on Z^^^{Tn). In the case ol K — ?, this equation looks 
like 



27 



+9 M-^"-^/^ < 



n>-2 



n>-2 



^(3/. - 2)r3,_2>v£U + E 



a+6=3Ti ^^3a+l 
a,6>0; n>-3 



36-3 



> + 



+9 M-^"-2/3 , 



n>-2 



a,+b=3n ^ 3a+l 
a,6>0; n>-3 



= 0, 



(4.36) 



(see (1.27), (1.28) for explicit expressions for W-operators) . This relation is valid for any 
value of N , and in this sense is an example of identity, which behaves smoothly in the 
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limit iV — > cxo. 

— If = cxD all the quantities Tre{M)M~P with given p but varying e(M) become 
independent, and (4.36) may be said to imply (1.26). 

In what follows we restrict ourselves to an illustrative calculation: we prove that (4.36) 
holds for A?" = 1 (i.e. when there are numbers instead of matrices). Even this calculation 
is long enough. Additional technical details, necessary to deal with matrices rather than 
ordinary numbers, are exhaustively discussed in the previous section 4.2. 

Example. When (4.32) with A' = 3, A^ = 1 is substituted into (4.33) and the rule 
(4.34) is used, we obtain {Z = 2'P^) : 




dZ 



a 



+ E 



^ ,,n+l 

„>1 P ^-^n At 

1 d^z 

L^„^„27/i'+™+"9T,9T™9T„ 
n2 + 12n + 39 1 dZ 21 



1^ 1 d'^Z 



n + m + 



^ n + 4 1 dZ 7 ■ 

1 d^Z 



+ 



18 dT^dT^ 



27 



/x"ar„ 27^ 



(4.37) 



Our purpose now is to compare this expression with the l.h.s. of (4.36). In order to 

understand the structure of (4.36) let us note, that the first item (with W^^^) in (4.36) 

1 d^Z 

contains the contribution — —777^-7:7^——- , iust the same as in (4.37), but only under 
restriction I -\- m -\- n = mod 3. However, no such restriction is imposed in (4.37), 

and in order to restore the equivalence between (4.36) and (4.37), one needs to add the 

terms with ^(yy^^) Z)/dT to the l.h.s. of (4.36), which add the missing contributions 

1 d^Z 

&Tm^&T^ "''^ / + m + n = 1 , 2 mod 3. 

If we analyze the terms without a in (4.37), it is important that a appears also in the 
shift (4.35): 



nf^^ = nTr, - 3a6n,4, (4.38) 



so that at the same time we substitute all f's by T's in (4.36), (1.26)-(1.28). Then 

(AT- i)(3a:- i)(5a:- i)...((2a:- i)a:- 1) ^ 28^ „x , 

the term ^ '-^ , '- '-Z = —Z in 4.37 should be com- 

(2A^)^ 27 ^ ^ 

pared to the contributions without Z-derivatives to (4.36). These come from the negative 
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harmonics of W"-^^ and W*^^-* operators. Namely, so modified wi^g contains (2T2)^ + 

3(Ti)^(4T4) = — , >V_3 contains (Ti)^ = — , and there is no constant term in Wq . As 
A* A* 

to Virasoro operators, there is ^ • 2(ri)(2T2) = ^ in W^] and ~ 24K ^ ^ ^ 
in W^^^ . Also the contribution of interest to dW^^l/dTi is ^ • 2(2r2) = , and 

19 11 

;(2) /OT-i • _ _ Arrp \ _ Ai /'OT-' _^-iy\;(2) v.;k,,+qo _ . Of^AfO^A"^ — 



to dyV)_i/dT2 is - • 4(Ti) = — . Also {2T2)W]:i contributes - • 2(Ti)(2T2) ^ , 



, (4T4)>V(i^ - (4T4) ■ - = — and {Ti)W^o - (Ti) ■ - = — . Putting this all together we 
obtain the coefficient 



— {[4 + l]+9[- + - + - + -] + 9[- + - + - + -]] + 0(a) = — + 0(a) 
27 V ^ ^3 3 9 9^ ^3 3 9 9^/ ^ ' 27 ^ ' 

in front of Z in (4.36), in accordance with (4.37). In order to reproduce the term 

^ ^ 2AK ^ ^ ~~ — ^^^"^-^ ~ g*^ ^ (4.37) it is necessary to restore the shift 

(4.38) of T4- variable in (4.36). There are also contributions with r| in W*^^^ and in 
(3A; — 2)T'3fe_2W3^^+3„ , which are responsible to occurrence of the a^-term in (4.37). 

This illustrative comparison of (4.36) and (4.37) is enough to recognize the proper 
structure of (4.36). After it is found out it is easy to ensure that all the remaining 
contributions to (4.36) and (4.37) also coincide. Of course, in order to check the trace- 
structure of (4.36) our simple example of = 1 is not enough — the analogue of detailed 
consideration of sect.4.2 is required. Such detailed derivation, as well as consideration of 
the case of X > 3 , is beyond the scope of this paper. 



5 Conclusion 

To conclude, we presented enough evidence that the Generalized Kontsevich Model (1.2) 
interpolates between all the multimatrix models, while preserving the property of inte- 
grability and the >C_i-constraint. This makes GKM a very appealing candidate for the 
role of a theory, which could unify all "stable" {i.e. with c < 1) bosonic string models. 
However, the study of GKM was at most originated here. Let us list a set of problems, 
which seem interesting for the future development. 
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1) Though integrable structure and C-i- constraint in GKM have been discussed more 
or less exhaustively, the situation with generic W-constraints and their implications re- 
mains less satisfactory. Of course, they may be derived from £_i-constraint, as suggested 
in refs.[|ll, but still it is desirable: 

— to complete the derivation of W-constraints from Ward identity (1-23) for = cxd, 
as suggested in sect. 3, for > 3 ; 

— to prove the equivalence of these W-constraints (which look like (1.26)) to conven- 
tional constraints (1.29); 

— to describe the subvariety in Grassmannian (at infinity of the Universal module 
space), specified by the universal C^^^- constraint; 

— to take double-scaling continuum limit of the properly reduced W-constraints in 
discrete multimatrix models (which look like Wg-p^"* [Mi] = W^^q^ [M2] for the 2-matrix 



case [W]) in order to derive the constraints (1.31) for yV]^'^^ (as it was done in ref. |T6 
for the 1-matrix case); 

— also the problems discussed in sect. 3. 4 should be added to this list. 

2) We tried to argue, that particular {K — l)-matrix model arises from GKM for a 

particular choice of potential: V(X) = Vk{X) = const ■ X^^^. Particular {CFT + 2d 

gravity)-model is specified by additional adjustment of M-matrix in such a way, that all 

Tn = 0,n^l, p + K (for c = 1 ). These particular points in the space 

pK 

of parameters {V{X), M} should be interpreted as critical points of GKM. The relevant 
questions are: 

— what is the proper criterium, distinguishing critical points in terms of GKM itself; 

— are there any other critical points? 

One can even hope, that the answer to the last question is positive and some other 
critical points may be identified with non-bosonic string models. Let us remind that the 
stability criterium like c < 1 (implying the absence of tachionic instabilities, leading to 
desintegration of Riemann surfaces through creation of holes), which in bosonic case im- 
plies the restriction d < 2 for the space time dimension, in the case of = 2 superstrings 
turns into a much more promising condition: d < 4 . 

3) A separate class of problems is related to identification of (double scaling limit) 
c = 1 models in the theory of GKM. This could lead to a more natural formulation of 
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GKM, in particular, help to restore the symmetry between arguments V{X) and M of 
Z^^J^ [M] . Keeping in mind the topological implications oi K = 2 model it is also natural 
to ask, 



what corresponds to Penner model 38, 39 



4) The last group of problems concerns the meaning and generalizations of GKM. First 
of all, both formulations in terms of topological (like ^ 41, 39 1) and conformal 

(in the spirit of is required in order to understand explicitely the connection 

to Liouville theory. Second, the hope that just GKM (1.2) is enough to include N = 
2 superstrings may be too optimistic, and then one needs to look for more universal 
models. Third, in any case one should search for a more invariant (algebro-geometric) 
interpretation of (1.2): GKM is already much simpler than original multi-matrix models 
(it is a 1-matrix model with a smooth continuum limit), but still it is not simple enough 
to be appealing as a fundamental theory. The main aesthetic problem is the assymetry 
between M and V{X) and, perhaps, the very relevance of matrix integrals. Fourth, 
GKM does describe interpolation between sufficiently many string models, but it is not 
a dynamical interpolation: one may still need to look for an "of-shell" version of GKM, 
involving a sort of integration over V(X). This problem can be after all related to the 
question about critical points. In particular, it may be interesting to find out something 



similar to the "off-shell" actions of refs.[17, 44, 45], which reproduce Virasoro constraints 
in the ordinary 1-matrix models as dynamical equations of motion, for the case of GKM. 



Possible links between GKM and double-loop algebras (in the spirit of |^) also deserve 
investigation. 
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